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ABSTRACT 


A  high-aspect-ratio  planing  surface  gliding  on  a 
stream  of  an  infinitely  deep,  incompressible,  inviscid  and 
gravity-free  fluid  is  treated.  This  complicated  problem  is 
decomposed  into  two  relatively  simpler  boundary- value 
problems . 

The  near-field  boundary-value  problem  is  valid  only  in 
the  neighborhood  of  the  planing  surface.  The  problem  is 
solved  by  the  classic  hodograph  method.  The  second-order 
inner  problem  is  also  shown  to  be  a  plane,  irrotational 
flow  and  the  solution  is  obtained  by  following  the  same 
procedure  as  given  in  the  first-order  inner  solution. 

The  far-field  boundary-value  problem  is  valid  only 
far  away  from  the  planing  surface.  The  first-order  outer 
solution  is  shown  to  be  a  trivial  uniform  flow.  The  outer 
velocity  potential  is  defined  in  the  whole  space  by  har¬ 
monic  continuation.  The  second-order  solution  is  then 
shown  to  be  similar  to  a  lifting  line  solution. 

The  unknown  strength  of  singularities  is  obtained  by 
matching  of  the  velocity  potential.  Then,  a  matching  of 
the  free-surface  deflection  provides  a  height  reference  for 
the  planing  surface.  The  location  of  che  planing  surface 
with  respect  to  the  undisturbed  free  surface  is  uniquely 
defined. 

In  order  to  obtain  a  unique  second-order  solution,  it 
is  necessary  to  solve  the  third-order  outer  solution.  The 
detail  of  this  solution  is  presented. 


ii 


A  numerical  solution  for  a  planing  plate  of  arbitrary 
angle  of  attack  is  presented.  A  downwash  correction  is 
also  included. 

It  is  shown  mathematically  that  the  present  theory 
can  be  applied  to  V-shape  or  general  shape  planing  surfaces 
with  curvature  in  the  spanwise  direction. 


CONTENTS 


Page 

ABSTPvACT .  ii 

SYMBOLS . . .  Vi 

CHAPTER 

I .  INTRODUCTION 

1.  Introduction  and  Review .  1 

2.  Scope  and  Nature  of  the  Solution .  5 

II.  FORMULATION  OF  THE  PROBLEM 

1.  Statement  of  the  Problem .  8 

2.  Boundary-Value  Problem  for  the 

Planing  Plate . 9 

3.  The  Far-Field  Boundary-Value  Problem...  15 

4.  The  Near-Field  Boundary-Value  Problem..  18 

III.  FIRST-ORDER  SOLUTION 

1.  First-Order  Inner  Solution .  24 

2.  First-Order  Outer  Solution .  33 

3.  First-Order  Matching .  36 

IV.  SECOND-ORDER  SOLUTION 

1.  Second-Order  Outer  Solution .  40 

2.  Intermediate  Step  Matching .  45 

3.  Second-Order  Inner  Solution .  49 

4.  Second-Order  Matching .  57 

V.  FREE-SURFACE  MATCHING 

1.  Introduction .  60 

2.  First-Order  Free-Surface  Matching .  61 

A.  Near-Field  Free-Surface  Solution...  61 

B.  Far-Field  Free-Surface  Solution....  62 

C.  First-Order  Free-Surface  Matching..  63 

3.  Second-Order  Outer  Free-Surface 

Solution .  64 

4.  Intermediate  Free-Surface  Matching .  65 

5.  Second-Order  Inner-Surface  Solution....  67 

6.  Second-Order  Free-Surface  Matching .  67 


IV 


CONTENTS  (CONT'D) 


Page 

VI.  HYDRODYNAMIC  LIFT  CALCULATION 

1.  Chord  Length  of  the  Planing  Surface....  69 

A.  First-Order  Chord  Length  Expression  69 

B.  Second-Order  Chord  Length 

Expression .  70 

2.  Wetted  Length  of  the  Planing-Surface. . .  70 

3.  Jet  Direction .  71 

4.  Lift  Force .  72 

A.  First-Order  Lift  Force .  73 

B.  Second-Order  Lift  Force .  74 

5.  First-Order  Numerical  Calculation .  75 

A.  Cusp  Shape  Planing  Plate .  77 

B.  Elliptical  Shape  Planing  Plate .  77 

6.  Second-Order  Numerical  Calculation .  86 

VII.  THIRD-ORDER  OUTER  SOLUTION 

1.  Introduction .  90 

2.  Third-Order  Outer  Boundary-Value 

Problem .  90 

3.  Third-Order  Outer  Solution .  91 

4.  Intermediate  Velocity  Potential 

Matching . 94 

5.  Third-Order  Free-Surface  Expression....  98 

VIII.  WEDGE  SHAPE  PLANING  SURFACE 

1.  Introduction .  101 

2.  Formulation  of  the  Problem .  101 

3.  Far-Field  Boundary-Value  Problem  and 

Solution .  103 

4.  Near-Field  uoundary-Value  Problem  and 

Solution .  103 

5.  General-Shape  Planing  Surface .  105 

IX.  CONCLUSION . 108 

REFERENCES .  Ill 

APPENDICES .  114 


V 


SYMBOLS 


Primed  variables  are  dimensional;  unprimed  variables  are 
dimensionless 


Aq^Ai 

AR 


^O'^l 


“o'^l 


H 


L 

£ 


Mq^Mi 


N 


No,Ni 

P 


Rq^Ri 


constants  in  inversion  formula 
aspect-ratio 

jet  thickness (outer  variables) 

=ao+ai 

jet  in  C"plane 

stagnation  point  in  ^-plane 


complex  potential 

height  of  the  trailing  edge  above  the 
undisturbed  free  surface  (inner) 

height  of  the  trailing  edge  above  the 
undisturbed  free  surface  (outer) 

local  chord  length  (inner) 

local  chord  length  (outer) 


constants 

unit  vector  in  a  plane  z=constant, normal 
to  body  contour  in  that  cross-section 

free-surface  elevation  (inner) 

pressure 

atmospheric  pressure 

constants 

free-surface 


vi 


SYMBOLS  (CONT'D) 


constant  unifoiro  velocity 

complex  velocity 

rectangular  coordinates  (outer) 

rectangular  coordinates  (inner) 

complex  variable 

angle  of  attack 

jet  angle 

strength  of  singularities 

jet  thickness  (inner) 

small  parameter 

complex  variable 

free-surface  elevation  (outer) 

mass  density 

lift  force  per-unit-span 

velocity  potential  (inner) 

velocity  potential  (outer) 

stream  function 

imaginary  unit 

real  part 


vii 


I.  INTRODUCTION 


I.  1.  Introduction  and  Review 

With  the  appearance  of  seaplanes  and  high  speed 
planing  boats,  the  problems  associated  with  planing  surfaces 
became  practical  problems.  It  is  interesting  to  note  that 
the  first  experiment  with  a  flat  planing  surface  was  carried 
out  in  1912,^  which  was  followed  by  a  series  of  both  theoreti¬ 
cal  and  experimental  investigations. 

At  rest  and  low  forward  speeds  the  planing  boat  is 
supported  by  buoyancy  force.  As  speed  increases,  the  center 
of  gravity  gradually  begins  to  rise.  The  center  of  pressure 
on  ti.  bottom  shifts  toward  the  nose  and  increases  the  trim 
by  the  stern.  With  further  increase  in  the  speed,  the  wetted 
area  decreases,  the  mean  draft  of  the  hull  becomes  so  slight 
that  most  of  its  weight  is  supported  by  the  dynamic  lift 

2 

acting  on  its  bottom,  and  the  hydrostatic  force  is  very  small. 

This  condition  is  called  ’planing  and  the  bottom  surface  which 

is  in  contact  with  water  is  called  the  planing  surface. 

3 

Wagner,  in  his  fundamental  paper,  showed  that  for  infini¬ 
tely  small  angles  of  attack,  with  irrotational  and  gravity- 
free  motion,  the  flow  on  the  lower  surface  of  a  planing 
surface  is  comparable  with  that  on  the  lower  surface  of  an 
airfoil.  The  lift  of  a  planing  surface  is  approximately  half 

as  great  as  that  of  an  infinitely  thin  airfoil  whose  plan 

4  5  6 

corresponds  to  the  wetted  surface.  Sottorf  '  and  Sambraus 
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indicated  a  satisfactory  agreement  for  moderately  small 
iiKiles  of  attack  between  their  test  data  on  flat  planing  sur¬ 
faces  and  Wagner's  theory.  However,  Wagner's  method  is  approx¬ 
imate  and,  strictly,  only  applies  when  the  angle  of  incidence 
of  the  gliding  plate  to  the  stream  is  moderately  small;  more¬ 
over,  his  solution  requires  that  the  direction  of  the  jet  or 
spray  which  is  formed  should  be  parallel  to  the  length  of 
the  plate.  In  the  middle  1930's,  A.  E.  Green  '  presented 
several  papers  on  this  subject.  He  treated  a  two  dimensional 
plate  gliding  on  the  surface  ,>f  a  gravity-free  stream  of 
infinite  depth.  He  was  able  to  obtain  a  complete  solution 
of  a  two  dimensional  gravity-free  planing-plate  problem, 
satisfying  the  non-linear  free-surface  boundary  conditions 
exactly. 

9 

However,  his  solution  is  not  unique.  There  is  another 
anomaly  in  Wagner's  and  Green's  theories  of  infinitely  deep 
water:  their  theories  imply  that  the  free-surface  elevation 
infinitely  far  upstream  and  downstream  is  depressed  an 
infinite  amount. 

The  principal  difficulty  of  a  free  surface  problem  is 
that  both  the  dynamic  and  kinematic  free-surface  boundary 
conditions  are  non-linear.  Worse  than  that, the  free 
surface  on  which  the  conditions  are  to  be  applied  is  not 
known  a  priori;  it  has  to  be  found  as  part  of  the  problem 
solution.  To  overcome  this  difficulty,  another  approach 
is  found  in  the  literature.  Assume  the  disturbance  is 
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small;  then  the  problan  can  be  linearized.  It  is  assumed  that 
the  error  is  of  high  order  if  the  free  surface  boundary 
conditions  and  the  body  condition  are  applied  on  the  undisturb¬ 
ed  free  surface  instead  of  on  actual  positions  of  free  surface 
and  body  surface.  Then  the  higher  order  terms  in  the  free 
surface  boundary  conditions  are  neglected  compared  to  the 
linear  term.  In  this  way  the  two  dimensional  boundary-value 
problem  is  reduced  to  an  integral  equation  which  was  first 

solved  by  Sretensky^®  later  by  Maruo^^  and  Cumberbatch. 

13 

Lamb  considered  the  two  dimensional  flow  due  to  the  appli¬ 
cation  of  a  pressure  distribution  on  the  surface  of  a  stream. 

He  computed  two  simple  pressvire  distributions  for  which  the 
integrals  can  be  evaluated  to  give  the  shape  of  the  planing 
surface.  Such  linearized  solutions  give  a  good  approximation 
to  the  expected  physical  behavior,  except  in  the  neighborhood 
of  the  plate,  where  the  flow  near  the  stagnation  point  and 
the  spray  sheet  can  not  be  regarded  as  small  perturbations 

of  the  uniform  streeim. 

14  15 

Rispin  and  Wu  used  the  singular  perturbation  method 
to  Folve  the  two  dimensional  non-linear  planing  problem  by 
including  gravity  effects.  They  obtained  a  complete  solution 
which  includes  an  inner  expansion  and  an  outer  expansion. 

Both  anomalies  in  Wagner's  and  Green's  classic  solutions 
were  removed. 

In  the  experimental  approach.  Shoemaker^®  made  a  test 
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for  V-shaped  as  well  as  flat-plate  planing  surfaces.  Korvin- 

Kroukovsky^^  gave  an  empirical  foimula  for  lift  calculation. 
13  19 

Shuford  *  gave  a  review  of  planing  theory  and  experiment 

of  rectangular  flat  plates.  Using  the  semi-empirical  formula 

the  effects  of  cross  section  and  plan  form  were  included. 

Progress  in  developing  a  three  dimensional  theoretical 

20 

approach  is  rather  slow.  Using  linearized  theory,  Maruo 

obtained  an  approximate  solution  from  the  integral  equation 

for  high-eind  low- aspect-ratio  planing-surface  problems. 

A  low-aspect-ratio  planing  surface  has  been  treated  as  a 

21 

slender  body  by  Tulin,  without  matched  expansions  and  by 
22 

Ogilvie,  with  matched  expansions. 

A  low-aspect-ratio  planing  surface  is  relatively  ineffi¬ 
cient  in  producing  lift  force;  more  area  has  to  be  provided 
for  a  given  lift,  with  a  resulting  penalty  in  lift-drag  ratio 

compared  to  that  of  a  high-aspect- ratio  surface  such  as  a 

23 

hydrofoil.  Clement  proposed  a  new  plan  ^orm  with  a  re¬ 
entrant  step  in  his  planing  boat  design.  Brown  and  Van 

24  ..... 

Dyck  made  an  experimental  investigation  on  this  high- 

aspect-ratio  re-entrant  step  planing  surface,  with  an  encoura 

.  25 

ging  result.  Current  development  of  practical  stabilizers 

is  proceeding  rapidly  enough  that  within  a  few  years  a  high- 

aspect-ratio  planing  boat  will  become  practical. 

In  his  experimental  studies  on  a  rectangular-plan-form, 

26 

high-aspect-ratio,  planing  plate,  Mottard  found  that  an 
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instability  quite  similar  to  the  flutter  of  an  airfoil  takes 
place  even  if  only  a  single  degree  of  freedom,  such  as  heave 
motion,  is  present.  This  is  a  rather  interesting  result. 

In  aerodynamics,  the  occurrence  of  flutter  normally  requires 

27 

that  two  vibrational  degrees  of  freedom  be  involved.  Ogilvie 
linearized  the  lifting-surface  problem  in  the  manner  of 
Wagner  and  solved  this  two  dimensional  hydrodynamic  problem 
in  the  presence  of  heave  motion.  Le  obtained  a  parameter  v, 
a  reduced  frequency,  which  characterized  the  occurrence  of 
instability.  His  theoretical  result  agrees  very  well  with 
Mottcird's  experimental  data.  In  order  to  solve  a  three 
dimensional  instability  problem,  a  three  dimensional  steady 
state  hydrodynamic  problem  has  to  be  solved  first.  This  is 
the  aim  of  the  present  v;ork« 

I.  2.  Scope  and  Nature  of  the  Solution 

The  trailing  edge  at  the  step  of  a  planing  boat  is  a 
straight  line.  The  contours  of  the  bottom  lines  of  the 
planing  surface  in  the  immediate  neighborhood  of  the  step 
closely  approximate  those  of  a  flat  plate.  As  the  speed 
of  the  planing  surface  increases,  the  wetted-aspect- ratio, 

AR,  increases  rapidly.  The  hydrostatic  force  and  friction 
force  become  negligible  compared  to  the  increased  hydro- 
dynamic  force.  In  the  present  study,  only  the  limiting 
high-speed  case  of  inviscid  flow  will  be  considered >  that 
is,  the  effects  of  gravity  and  viscosity  are  neglected. 
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In  order  to  apply  the  hodograph  method,  a  planing  surface 
without  camber  in  the  longitudinal  direction  is  considered. 

Let  the  aspect- ratio  AR  be  defined  as 

2 

ar  =  _  Span _  ^  _  Span _ 

Wetted  Planing  Area  Mean  Wetted  Chord  Length 

For  a  high-aspect  ratio  planing  surface,  it  is  implied 
that  the  maximum  dimension  in  the  z  direction  is  much 
greater  than  the  m£ucimum  dimension  in  the  cross-sections. 
Then  a  small  parameter  c  can  be  defined  as  the  inverse  of 
high-aspect  ratio: 

1  Mean  Wetted  Chord  Lengtii 

AR  ■  Span 

As  aspect- ratio  AR  becomes  larger  and  larger,  the  span 
grows  bigger  and  bigger  if  the  mean  wetted  length  is  fixed. 
At  its  limit  as  AR  approaches  infinity,  the  span  approaches 
infinity.  For  an  infinite  span  the  flow  will  be  two 
dimensional.  This  is  exactly  the  lowest-order  near-field 
description.  The  characteristic  length  will  be  the  chord 
length  at  each  section.  On  the  other  hand  if  e  becomes 
smaller  and  smaller,  the  scale  of  mean  wetted  length  becomes 
smaller  and  smaller  if  the  span  is  fixed.  At  its  limit 
as  e  approaches  zero,,  the  planing  surface  shrinks  to  a  line 
and  uniform  flow  is  recovered.  This  is  exactly  the  low¬ 
est-order  far-field  description.  The  characteristic 
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length  will  be  the  span.  There  are  two  different  reference 
length  scales  working  at  the  same  time.  It  suggests  that  a 
singular  perturtiation  method  can  be  applied  to  solve  this 
problon. 

In  the  present  %irork,  the  method  of  matched  asymptotic 

expansions  is  applied  to  solve  a  high-aspect- ratio  planing 

surface  probl^.  The  full  boundary-value  problem  is 

decomposed  into  near-field  and  far-field  boundary-value 

problems.  The  first-order  near-field  problem  with  respect 

to  the  small  parameter  c  is  found  to  be  a  plane  irrotational- 

flow  problem.  The  solution  is  obtained  by  the  hodograph 

method.  The  first-order  outer  solution  represents  a 

uniform  stream.  The  second-order  inner  problem  is  shown 

to  be  a  plane  irrotational-flow  problem  again.  The 

second-order  outer  problem  is  shown  to  be  similar  to  the 

28 

lifting  line  problem.  Following  the  matching  principle 
as  outlined  by  Van  Dyke,  unknown  parameters  in  the  near- 
and  far-field  problem  are  determined.  Finally,  a  unique 
solution  is  obtained  through  matching  of  the  free-surface 
deformation.  The  depth  of  submergence  can  thus  be  properly 
defined.  The  anomaly  in  the  Wagner's  and  Green's  gravity- 
free  infinite  depth  solution  is  thus  removed  by  the  present 
work.  Three  dimensional  effects  are  shown  to  provide  a 
height  reference.  Thus  present  theory  is  applicable  for 
heavy  loading  as  well  as  light  loading  on  the  planing  surface. 


II.  FORMULATION  OF  THE  PROBLEM 

II.  1.  Statement  of  the  Problem 

The  problem  treated  here  is  a  three  dimensional  flow 
generated  by  a  high-aspect-ratio  planing  surface  gliding  with 
a  constant  speed  U*  at  an  arbitrary  angle  on  the  free  surface 
of  an  infinitely  deep  weightless  fluid.  The  fluid  is  assumed 
to  be  inccmipressible  and  inviscid  and  the  flow  irrotational . 
Instead  of  a  planing  surface  moving  on  calm  water,  a  uniform 
stream  of  speed  U'  at  infinity  in  the  direction  of  positive 
x-axis  is  superimposed.  Then  a  steady  motion  is  obtained. 
Take  the  coordinate  system  with  the  origin  on  the  undisturbed 
free  surface  and  choose  the  x-axis  in  the  direction  of  the 
uniform  stream.  The  positive  y-axis  is  taken  vertically 
upwards  and  the  z-axis  is  perpendicular  to  the  xy-plane  as 
shown  in  Figure  2.1. 

Let  the  span  be  b*  and  the  planing  surface  located 

in  the  segment  |z|  £  b'/2.  Let  the  chord  length  of  the 

planing  surface  be  ll'(z').  The  function  il'(z')  must  be 

continuous  and  have  a  continuous  derivative,  at  the  tip  as 

elsewhere;  otherwise  singularities  will  occur  and  the 

28 

region  near  the  tips  will  need  a  special  treatment. 

Let  the  angle  of  incidence  be  a,  the  aspect  ratio  of  the 
planing  surface  be  AR,  and  the  small  parameter  be  e,  which 
is  the  inverse  of  the  aspect  ratio.  Let  the  vertical 
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distance  between  the  trailing  edge  and  the  undisturbed 
free-surface  level  be  h'(2').  To  simplify  the  analysis, 
it  is  assumed  here  that  the  trailing  edge  is  a  straight  hori¬ 
zontal  line.  This  assumption  implies  that  h'  is  independent 
of  z.  Therefore  the  uncambered  planing  surface  beccmies  a 
planing  plate.  In  Chapter  8,  this  restriction  will  be 
relaxed  so  that  the  bott<xn  of  the  planing  surface  can  be 
V-shaped  or  even  curved  in  the  spanwise  direction. 

Assume  that  the  angle  of  incidence  a,  the  depth  of 
sutoergence  h*  2md  the  shape  of  the  chord  length  distri¬ 
bution  V  (z*)  are  given.  The  solution  of  the  hydrodynamic 
problem  of  the  planing  surface  is  to  be  determined. 

II.  2.  Boundary  Value  Problem  for  the  Planing  Plate 

The  fluid  is  assumed  to  be  ideal  and  the  flow  irrota- 
tional.  There  exists  a  velocity  potential  <!»' (x' ,y' ,z') 
which  satisfies  the  Laplace  equation 

(L)  =0  in  the  fluid  region 

Assume  that  the  free  surface  deflection  can  be  specified 
by  the  equation 


on  the  free  surface 


y'  =  n'  (x' ,2' ) 
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The  function  n‘  is  unknown;  it  is  a  part  of  the  problem 
to  be  solved. 

Tliere  are  two  conditions  to  be  satisfied  on  the  free 
surface.  Let  the  velocity  be  the  positive  gradient  of 
the  velocity  potential,  q'  =  V<p* ,  where  q'  is  the  velocity 
vector.  The  first  condition  is  called  the  dynamic 
boundary  condition.  It  states  that  the  pressure  is  constant 
on  the  free  surface,  equal  to  the  atmospheric  pressure: 


p'  (x* ,n' (x* =  p^’ 


on  Sj 


The  Bernoulli  equation  gives 


P*  .  J,  4. 


p;  1  2 


where  p^  is  the  atmospheric  pressure,  p'  is  the  local 
pressure  xn  the  fluid,  p'  is  the  density  of  the  fluid.  An 
assumption  is  made  that  no  disturbance  exists  far  upstream 
at  infinity.  On  the  free  surface,  the  Bernoulli  equation 
becomes : 


(FI) 


=  u' 


on  Sg 


The  second  condition  on  the  free  surface  is  called 
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the  kinematic  boundary  condition.  It  states  that  a 
particle  once  in  the  free  surface  remains  in  the  free 
surface.  Mathematically,  it  can  be  stated: 


Dt 


^  I.r,'(x',z')  -  y'l  =  0 


on  S 


f  , 


where  ^  denotes  the  substcuntial  derivative.  The  above 
equation  gives: 


(F2)  -  <|)y,  +  =0  on 


Equations  (FI)  and  (F2)  are  two  non-linear  conditions  to 
be  satisfied  on  an  unknown  boundary 

Next,  a  body  boundary  condition  must  be  satisfied. 

A  body  surface  might  be  defined  by  the  equation 

s'  =s'(x',z')  +  y'  =0 

The  kinematic  boundary  condition  on  the  body  would  be 


4’'  |S'  ,  - 
x'  x' 


d)' ,  +  4>'  ,s'  =  0 

y  z  z' 


on  y'  =  -  s' (x' ,z') . 


For  the  planing  plate,  for  which; 


i 

\ 

I 

f 
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s'  =  X '  tan  a  +  y '  =  0 , 

tne  body  condition  is  : 

(B)  4)J_,  sin  a  +  <p^,  cos  a  =  0  on  the  plate 

In  airfoil  theory,  a  Kutta  condition  is  applied  to 
give  a  finite  velocity  at  the  trailing  edge.  In  the  planing 
problem,  this  condition  requires  that  the  flow  separates 
smoothly  from  the  the  trailing  edge  of  the  planing  surface. 

Finally,  there  is  a  radiation  condition  to  be  satisfied 
to  guarantee  a  unique  solution.  This  condition  states  that 
there  is  no  disturbance  far  upstream. 

The  planing  plate  problem  is  thus  mathematically 
formulated.  The  governing  equation  is  the  Laplace  equation 
subject  to  two  free-surface  boundary  conditions,  a  body 
condition  and  a  radiation  condition. 

It  is  more  convenient  to  work  on  the  problem  by  making 
the  variables  and  the  equations  dimensionless. 

Let  the  span  be  normalized  to  be  of  length  2  and 
located  in  the  segment  |zj  ^1.  And  define  the  following 
variables  : 

I  b'  _  o 

^  ■  bV2  ^ 

x'  ,  y'  ,  z' 

X,y,Z  -  ^Tj2 
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q  = 


(u  V  w  =  V  •  w* . 

y,  vu,v,w  u*  »  u»  »  u*' 


n’  h' 

■  bV2  '  ^  "  bV2  '  ^  "  bV2 


p  =  — ^ 
p'U'^ 


.  _  4)' 

^  ■  “tPbV2 


The  full  boundary-value  problem  in  dimensionless  form 
becomes : 


(L) 

(FI) 

(F2) 


(B) 

(R) 


(b  +  (t»  +  <t>  =0 

^xx  ^yy  ^zz 


2  2  2 
(|)  +  <t>  ^  +  (|>  =1 

X  ^y  z 


♦x"*  -  ♦y  ♦  ♦z''x  =  “ 


<|)j^sin  a  +  4iyCOS  a  =  0 


<l>x  =  1  ;  <l>y  =  <l>2  =  0 


in  the  fluid  region 


on  Sj 


on  Sj 


on  plate 


as  x-»— y-*-a> 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 


Now,  consider  the  flow 
plate.  There  is  an  incident 
uniform  in  the  x-direction. 
velocity  potential  is  x.  As 
normal  to  the  incident  flow. 


field  far  away  from  the  planing 
flovj  which,  at  infinity,  is 
In  the  absence  of  the  body ,  the 
£-*■0,  the  body  shrinks  to  a  line 
so  the  first  term  in  the  far 
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field  description  represents  simply  a  uniform  stream, upon 
which  a  line  of  higher  order  singularities  is  superimposed. 

On  the  other  hand,  consider  the  flow  field  in  the  neigh¬ 
borhood  of  the  planing  plate.  The  flow  in  the  whole  vicinity 
of  the  pi  ining  surface  can  not  be  regarded  as  a  small 
perturbation  of  the  uniform  stream  for  arbitrary  angle  a 

As  aspect-ratio  AR  approaches  infinity,  the  flow  in  the 
near  field  acts  as  if  £  is  fixed  and  the  scale  of  span  b 
approaches  infinity  at  each  section.  In  this  sense,  it  is 
similar  to  the  flow  of  infinite  span  and  the  flow  field 
is  two  dimensional  to  the  lowest  order.  So  the  singular 
perturbation  method  can  be  applied  to  reduce  the  original 
complicated  free-surface  problem  into  two  relatively 
simpler  problems.  The  outer  solution  is  valid  only  far 
away  from  the  planing  surface,  while  the  inner  solution 
is  valid  only  near  the  planing  surface.  An  approximate 

h 

solution,  valid  everywhere,  in  the  form  of  a  power  series 

in  e,  will  be  obtained  by  calculating  successive  solutions 

of  the  outer  and  inner  problems  and  matching  them  according 

29 

to  the  matching  principle: 

m-term  inner  expansion  cf  (n-term  outer  expansion) 

=  n-term  outer  expansion  of  (m-term  inner  expansion).  (2.6) 


*  The  power  series  contains  "logarithmic"  terms.  However, 
log  e  will  be  treated  as  0  (1)  as  suggested  by  Van  Dyke.^^ 
Otherwise  the  Van  Dyke  rule  leads  to  contradictory  results. 
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Here  in  and  n  are  any  two  integers.  This  matching  is 
necessary  in  order  to  obtain  inner  boundary  conditions 
for  the  outer  solution  and  vice  versa. 


II.  3.  The  Far  Field  Boundary- Value  Problem 

Let  X,  y,  z  be  the  outer  variables  which  correspond 
to  the  natural  variables.  Let  ((>  be  the  velocity  potential 
for  the  outer  problem.  The  xz-plane  is  located  on  the 
undisturbed  free  surface  with  origin  located  at  the 
intersection  of  the  planing  plate  and  the  undisturbed  free 
surface,  as  shown  in  Figure  2.2.  The  outer  problem 
corresponds  to  the  flow  field  far  away  from  the  plate. 
There  the  detail  of  the  planing  plate  and  the  Kutta 
condition  become  meaningless,  while  the  Laplace  equation, 
dynaii'ic  and  kinematic  free-surface  conditions,  and 
radiation  condition  must  still  be  satisfied.  The  boundary 
value  problem  for  the  far  field  can  thus  be  formulated  as 
following : 


(L)  (j)  +  (j)  +  (j)  =  0  in  fluid 

'  '  ^xx  ^yy  ^zz 

(FI)  +  <()2^  =  1  on  Sj 

(F2)  (}>^  Tljj  ~  <l>y  +  <f>z  =  0  on  Sj 


(2.7) 

(2.8) 

(2.9) 


(R)  ^  ^  '  ‘^y  ^  ^  ^  x-^-®,  y^-»  (2.10) 
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Next,  assvune  that  there  exists  an  asymptotic  expansion 
valid  in  the  far  field: 

H 

♦  (XfYfZ)  r  4'_(x,y,z;e) 

n=0  ^  (2.11) 

where 

as  e-*-0  for  fixed  x,y,2 

Also  assume  the  existence  of  a  complementary  asymptotic 
expansion  for  the  free-surface  elevation 

N 

n  '''  I  n_(x,2;e)  (2.12) 

n=0  " 

where  *  o  (n^^) 

as  e-^-O  for  fixed  x,2 

Substitute  equations  (2.11)  and  (2.12)  into  equations 
(2.7)  -  (2.10)  to  give; 

*^0xx  ^  *^0yy  ^  *^0zz  ^  *^lxx  "^lyy  *^lzz 

‘*’2xx  ‘^’2yy  ^2zz  +  ••••=  0  (2.13) 

♦ox^  ^  ♦oy^  *  ♦oz'  *  ♦oj,  ♦ij, 

+  2'»oz*lz  +  ♦ix^  *  ^  ♦ox  ♦2x  +  ♦ly^+  ^♦oy*2y 
+  *12^  "  ^^Oz+lz*  2''’1x+2x  ^♦ly^2y  ^‘'’lz*2z 

+......=  1  on  Sj  (2.14) 
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(F2)  '^ox^Ox  ■  ^Oy  ^Oz^Oz  ^Ox^lx  ■*■  '^’ix'^Ox  “  *^ly 

^  '^Oz^lz  ^  ♦iz'^Oz  ^  hx^lx  ^  ♦2x''0x  “^Ox  ^2x 

-  ‘**2y  ^  '^Oz  '^2z  ^  ^Iz  ^Iz  ^  ♦2z  ^Oz  ^  ••*'0 

on  S^  (2.15) 

1 

^  as  x-*^-"®,  y-^-«>  (2.16) 
0 

Equations  (2.13)  -  (2.16)  represent  the  mathematical 
expression  for  the  outer  boundary- value  problem. 

II.  4.  The  Near  Field  Boundary-Value  Problem 

In  the  neighborhood  of  a  high-aspect  ratio  planing 
plate,  the  flow  acts  as  if  it  is  nearly  two  dimensional. 

In  order  to  investigate  details  of  this  flow  field,  one 
thus  shall  have  to  become  more  and  more  nearsighted  as 
the  limit  is  approached.  Mathematically  the  x-and  y- 
coordinates  have  to  be  magnified  in  order  to  have  flow 
variables  change  on  a  reasonable  scale.  Let  x,  y,  z 


(R)  <tox  ♦lx  ♦2x 


♦oy  ♦ly  ♦2y  ^ 


♦oz  ♦iz  ♦2z 
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be  the  natural  variables. 

Let  X,  Y,  Z  be  the  magnified  inner  variables  with 
origin  at  the  trailing  edge,  *  tlie  inner  velocity  potential. 
Then 


X  = 


X  -t-  h  cot 
e 


a 

f 


Z  =  z  (2.17a) 


as  shown  in  Figure  2.3 


J 


Also  define 


H 


(2.17b) 


A  physical  quantity,  such  as  velocity,  should  not  be  changed 
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by  magnified  variables: 


llie  inner  problem  is  valid  only  in  the  neighborhood 
of  the  planing  plate.  The  radiation  condition  is  lost. 
However,  the  Laplace  equation,  dynamic  and  kinematic  free- 
surface  conditions,  body  boundary  condition  and  the  Kutta 
condition  must  still  be  true.  The  boundary -value  problem 
for  the  near  field  can  thus  be  formulated  as  follows: 


(L)  +  e  «2Z  “  ^  fluid  (2.18) 


(PI)  ^ 


on  S, 


(2.19) 


(F2)  ±  t  =  0 


(2.20) 


(B)  sin  a  +  "ty  cos  a  =  0  on  body  (2.21) 


(K)  The  flow  separates  smoothly  from  the  trailing 
edge 


Next  assume  the  existence  of  the  following  asymptotic 
series  with  respect  to  e 
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♦  '^n=0  ♦n<X,Y,Z;e) 

where  ♦  =  o{*  ) 

n+1  n 


(2.22a) 


e  0  for  fixed  X,  Y,  Z 


N  -v  ZN  (X,Z;e) 
n-0 


where  N  =  o(N  ) 
n+i  n 


e  0  for  fixed  X,  Z 


(2.22b) 


Substitute  equations  (2.22a)  and  (2.22b)  into  equations 
(2.18)  -  (2.21)  to  give: 


'^OXX  IXX  2XX  3XX 


■‘■'^OYY  “^lYY  '^2YY  "^lYY  ** 


(2.23) 
=  0  in  fluid 


^ox  ^%x^ix  ^"^ox^zx  "^IX  +  •  •  •  • 

2  2 

^OY  ^^OY^IY  "^lY  ^'^0y'^2Y  *  “  * 

^^OZ  ^^OZ^IZ  ^IZ  '^2Z'*’*‘“^ 


on  S, 


(2.24) 


e  ^'^OX^^OX  ■*■  ^OX^IX  "^IX^OX  ^IX^IX 

^0X**2X  ^  ••••)  "  ^'^OY  ^lY  ^2Y  “^SY  '*’*  * 


^^^0Z*^0Z  ^OZ^IZ  ^  '^17^07.  ■*■  ^17.**1 


1Z“0Z 


1Z”1Z 
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*  '0Z^2Z  * - *  ~  ° 


(2.25) 


(B)  +  'l>2x  ■*■  ••••)sin  a  + 

+  ■•■  ^2Y  *  ••••^cos  a  =  0  on  body  (2.26) 


This  is  the  non-linear  near-field  surface  problem. 

Now,  consider  the  lowest  order  of  this  inner  problem. 
According  to  the  definition  of  asymptotic  expansion, 
is  of  higher  order  than  equation  (2.23)  and  can  be 

discarded.  Therefore  the  possible  lowest-order  Laplace 
equation  is : 


"^OXX  ^OYY  ^  ^OZZ  ^ 


However,  according  to  the  definition,  <1*022 
scime  order  as  ^nvv*  Finally,  it  gives: 


(L) 


(b  +  =0 

OXX  OYY 


in  fluid  (2.27) 


Similar  argument  leads  to  the  following  boundary  condi¬ 
tions  . 


■'OX^  *  *0Y^  =  1 


on  Sj 


(2.28) 


(F2)  “^OX^OX  ”  ^^OY  ® 


on  S, 


(2.29) 
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(B)  sin  a  +  <!»gy  cos  a  =  0  on  body  (2.30) 

This  is  the  iowest-order  inner  boundary  value  problem, 
which  will  be  called  the  fipst- order  inner  problem  throughout 
the  entire  work. 


III.  FIRST-ORDER  SOLUTION 


III.  1.  First-Order  Inner  Solution: 

The  inner  boundary- value  problem  was  mathematically 
formulated  in  the  previous  chapter  and  then  was  expanded  in 
an  asymptotic  series  with  respect  to  e.  To  the  first 
order  the  boundary  value  problem  is: 


(L) 


*0XX  ‘^OYY  ® 


in  fluid  (3.1) 


(Fl)  ^  ^ 


on  S. 


(3.2) 


^OX^OX  “  ^^OY  ~  ° 


on  S, 


(3.3) 


(B)  sin  ot  +  <J>QY  cos  a  =  0  on  body  (3.4) 


(K)  Kutta  Condition 


The  governing  differential  equation  is  a  two  dimen¬ 
sional  Laplace  equation  at  each  section.  The  boundary 
conditions  are  all  two  dimensional.  Therefore,  the  first 
order  inner  problem  concerns  a  plane  irrotational  flow. 
The  dynamic  boundary  condition  (Fl)  states  that  the 
magnitude  of  the  velocity  on  the  free-surface  is  constant 
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and  known.  The  body  condition  (B)  states  that  the  direction 
of  the  flow  on  the  body  surface  is  known.  This  is  a  free 
streamline  problem  and  can  be  solved  by  the  classical  hodo- 
graph  method. 

It  is  noted  that  Z  coordinate  appears  only  implicitly 
in  this  first-order  inner  problem.  It  will  be  shown  in  the 
next  chapter  that  the  Z  coordinate  also  appears  only 
implicitly  in  the  second  order  inner  problem.  Therefore, 
in  order  to  stay  with  the  conventional  symbol,  Z  will  be 
defined  as  a  complex  variable  Z=X+iY  through  the  whole 
analysis. 

Experiments^^  indicate  that  near  the  planing  surface 
the  flow  shown  by  tests  presents  a  true  jet  or  spray 
character;  i.e.,  at  some  distance  ahead  of  the  planing 
body  the  water  surface  is  practically  undisturbed,  while 
immediately  forward  of  the  planing  body  the  water  is  thrown 
off  in  a  spray.  In  the  region  behind  the  planing  body, 
there  is  no  spray  separating  at  the  trailing  edge  at 
large  planing  speeds;  the  water  flows  off  smoothly  from 
the  trailing  edge. 

From  this  experimental  observation,  the  first-order 
inner  problem  can  be  stated  as  following:  An  infinitely 

long  flat  plate  is  held  at  an  arbitrary  angle  a  on  the 
free  surface  of  an  otherwise  uniform  stream  of  an  infinitely 
deep  weightless  fluid.  The  plate  is  of  length  L.  The 
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origin  is  located  at  the  trailing  edge  B  which  is  at  a 
depth  H  from  the  undisturbed  free  surface.  Y-axis  is 
vertically  upwards,  X-axis  is  in  the  direction  of  uniform 
stream.  The  plate  is  AB.  The  flow  ccanes  from  upstream 
infinity  and  divides  along  the  stagnation  streamline  which 
branches  off  at  the  stagnation  point  C.  Above  this 
streamline,  the  flow  shoots  off  as  a  jet  J.  Below  this 
streamline,  the  flow  separates  smoothly  from  the  trailing 
edge  B  as  shown  in  Figure  2.3.  From  equation  (3.1)  the 
velocity  potential  <*>Q(X,y)  is  harmonic  and  its  conjugate 
'i'Q(X,Y)  is  the  stream  function.  Define  the  complex 
potential  Fq(Z); 

Fg(Z)  =  <l'g(X,Y)  +  i'i*Q(X,Y) 

Let  Wq(Z)  be  the  derivative  of  Fq(Z)  with  respect 

to  Z; 


Wo(Z)  = 


dFQ(Z; 

dZ 


=  -  iVr 


Ox 


OY 


where  and  are  velocity  components  in  the  X  and  Y 
directions,  respectively.  Wq(Z)  is  called  the  complex 
velocity.  Both  Fq(Z)  and  Wq(Z)  are  analytic  functions. 
The  dynamic  free-surface  condition  (FI)  gives; 
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2  2 
"o  +  1 


on  Sj 


(3.5) 


The  kinematic  condition  gives: 


1  ^0  *0Y  '-•q 

^  <“■  ■  ♦ox  ■  “o 


on  S. 


(3.6) 


The  body  condition  gives: 


-tan  a 


(3.7) 


The  first-order  inner  problem  is  thus  expressed  in 
natural  variables.  This  problem  will  be  transformed  onto 
the  lower  half  of  the  auxiliary  ^-plane  as  shown  in  Figure 
3.2.  The  problem  is  then  solved  in  the  C“plane  and  finally 
the  solution  is  transformed  back  to  the  physical  plane  by 
the  inversion  formula. 

Let  the  potential  function  be  zero  at  C  and  the 
dividing  streamline  be  '1'q=0.  Consider  the  F^-plane  shown 
in  Figure  3.1.  IJ  is  the  upstream  free  surface.  If 
Uq  is  the  breadth  of  the  jet  at  a  great  distance  expressed 
in  the  outer  variables,  then  the  breadth  in  inner  variables  is 
.  Along  IJ,  the  stream  function  Tq  is 
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The  polygon  in  the  F^-plane  can  now  be  transformed 

onto  the  lower  half  of  C-plane.  Let  A,  B  correspond  to 

C  =  -1,  +1  respectively.  The  chord  length  is  represented 
by  the  segment  jcl  <  1  on  the  real  axis  n  =  0,  where 
C  =  r;  +  in.  Let  J,  C  then  be  mapped  onto  C  = 

The  free  surface  is  thus  mapped  onto  the  real  axis  too. 

The  diagram  of  the  C-plane  is  shown  in  Figure  3.1.  The 
Schwarz-Christof fel  transformation  gives: 


^ 

?  +  b. 


Integration  gives: 


-  <=0  -  ‘Vo' 
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This  is  Green's  solution.  (See  Milne-Thomson  ) 

Next,  define  the  mapping  function  Hq(c) 


Ho(rJ 


“"0 


TTb^+^  c+bo 


Then 


(3.9) 


~SZ~  ^O^^^dZ 


(3.10) 


From  equation  (3.10),  integration  gives 
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(3.11) 


The  origin  is  located  at  the  trailing  edge  to  give  the 
lower  limit  of  the  integral  to  be  1.  The  equation  (3.11) 
relates  the  fplane  and  the  physical  plane.  The  inner 
solution  found  in  the  C'plane  will  be  transformed  hack 
to  the  physical  plane  by  the  equation  (3.11) ,  which  is 
called  the  inversion  formula. 

The  equations  (3.5),  (3.6),  and  (3.7)  state  that  the 
flow  is  bounded  by  a  straight  rigid  boundary  (planing 
surface)  and  a  curvilinear  free  streamline  (free  surface) . 
The  complete  boundary  can  be  transformed  into  a  simple 
closed  polygon  by  a  Kirchhoff  function. 

Define  the  complex  variable  Qq  to  be 


Qq  =  -  log(WQ)  = 


{3.12a) 


Let  q  be  the  speed  in  the  flcft/  field.  Then 
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On  the  free  surface,  q=l,  which  gives: 


I  -  e-2’0  =  0 


on  S, 


(3.13a) 


On  the  body  surface. 


Gq  =  -  arg  Wq  =  -a 


Gq  =  TT  -  a 


on  CB 
on  AC 


(3.13b) 


Next,  express  these  boundary  conditions  on  the  real 
axis,  n=0,  of  the  ^-plane.  Equation  (3.13)  beccxnes: 


1  -  . 


O  (f)  =  {  ““ 

^  TT  -  a 


(n=o,U|>i) 


(n«0,CQ<c<l) 

(n=0,-l<5<CQ) 


(3.14) 


It  is  noted  that  the  kinematic  free-surface  condition  is 
automatically  satisfied  where  the  free  surface  on  the 
physical  plane  is  mapped  onto  the  real  axis  in  the 
segments  |cl  >  1  of  the  C”plane.  Equation  (3.14)  gives: 


To(5)  =  0 


So 


w. 


=  1 


(n=o,  K|>l) 


(n=o,  f;|>i) 


(3.15a) 
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And 


-a 


At  the  stagnation  point 


(n=0,Co<C<l) 

(.=0,-l<£<c„) 


"0<  =0 


)  =  0 


(3.15c) 


The  W^-plane  is  shown  in  Figure  3.2.  The  entire  fluid  region 
is  mapped  by  Wq(Z(5))  into  the  region  bounded  by  a  unit  half- 
circle  and  its  diameter.  The  diameter  is  the  image  of  the 
planing  surface,  on  which  the  direction  of  the  velocity 
vector  is  known,  and  the  semi-circle  is  the  image  of  the 
entire  free  surface,  on  which  the  magnitude  of  the  velocity 
vector  is  known.  Now,  map  the  W^-plane  onto  the  lower  half 
of  the  ?-plane  by  either  of  the  following  equivalent  forms: 


WQ(Z(t)) 


ia  ^"^0 

e  -  ' 

l-c^C+i/d-Cjj^)  (C^-1) 

e  '  ■  . .  . . 


(3.16) 


2  2  ^  /  r 

/l-^  is  defined  as  the  branch  of  (1-^  )  with  its 

branch  cut  on  the  real  c-axis  between  c=+l,-l. 
Equations  (3.9),  (3.11)  and  (3.16)  give: 
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Z(^)  = 


C+t>r 

71  (bQ+Cgi 

J^l- 


0  ■  '*  *  *"-0 
-  i/b.^-1  A-c^  log 


A-c^  loq(C+«^?^ 


log(C+A  -1) 

-/(bQ^^-l)  (c^-1) 


C+b, 


) 


(3.17) 


This  is  the  inversion  formula  which  gives  Z  as  a  function 
of  C. 

There  are  three  parameters  a^j  ,  b^  ,  and  c^  at  each 
section,  none  of  which  has  been  determined  yet  aside  from 
the  facts  that 

IcqI  <  1  ;  bQ  >  1 

Equations  (3.8)  and  (3.17)  provide  a  complete  first- 
order  inner  solution.  The  three  unknown  parameters  will 
be  determined  from  matching  to  the  far  field  problem. 

III.  2.  First-Order  Outer  Solution 

Consider  the  flow  field  at  a  great  distance  from  the 
planing  surface,  mathematically,  a  distance  of  order  unity 
with  respect  to  e.  The  boundary- value  problem  has  been 
formulated  in  section  II.  3.  With  the  same  kind  of  argument 
as  given  in  section  II. 4,  the  first  approximation  for  the 
outer  problem  is 
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The  free-surface  boundcury  conditions  (FI)  and  (F2) 
are  applied  on  which  is  unknovm  and  is  a  part  of  the 
problem  to  be  solved.  To  avoid  this  difficulty,  it  is 
assumed  that  Taylor  series  expansions  exist  for  at 
the  neighborhood  of  Hq-  From  (2.12),  is 

y  =  n(x,z)  =  Hq  +  n-  +  n2+  •••• 

It  is  then  assumed,  for  example,  that: 


4)Q^(x,n  (x,z)  ,z)  =  <I>Q^(x,nQ,z)  +  (Ti2+n2+  •  •  •)<l>oxy^*'’^0'^^ 


(ni+ri2+*  •  • ) 

+  2  '^Oxyy^^'’^0'^^  - 


However,  it  will  be  shown  in  the  later  part  of  this 
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section  that  the  first  term  of  the  outer  solution  represents 
a  uniform  stream.  It  gives  the  value  of  Hq  to  be  zero. 
Therefore,  it  is  assumed  here  that  nQ=0.  To  the  first- 
order  the  free- surface  conditions  can  be  applied  on  the 
undisturbed  free  surface  y^o  to  give: 


(FI) 

on  y-0 

(3.22) 

(F2) 

♦oy  =  “ 

on  y=0 

(3.23) 

For  a  high-aspect  ratio  planing  surface,  the  outer 
problem  corresponds  to  a  planing  surface  of  fixed  span 
and  vanishingly  small  chord  length.  As  AR  -►  the  planing 
surface  shrinks  to  a  line.  There  is  no  disturbance  far 
upstream.  Therefore,  it  is  plausible  to  assume  that  the 
outer  problem  can  be  linearized.  This  is  confirmed  in 
experiments^^  which  show  that  at  a  small  distance  ahead 
of  the  plcuiing  body  the  water  surface  is  practically 
undisturbed. 

This  physical  argument  shows  that  the  first  order 
solution  is  a  uniform  stream. 

(|>q(x,  y,  2)  =  X  (3.24) 

Mathematically,  it  can  be  shown  that  the  equation  (3.24) 
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satisfies  all  the  first-order  equations,  (3.18)  -  (3.23) 
trivially. 

III.  3.  First-Order  Matching 

In  the  inner  problem,  the  radiation  condition  is  not 
considered.  On  the  other  hand,  the  body  condition  is  not 
considered  for  the  far  field  solution.  These  missing 
conditions  can  be  recovered  through  the  matching  process, 
leading  to  a  unified  solution  valid  in  the  whole  fluid 
region. 

Let  m=l  and  n=l  in  the  equation  (2.6); 

1-term  inner  expansion  of  (1-term  outer  expansion) 

=  1-term  outer  expansion  of  (1-term  inner  expansion) 

(3.25) 

Each  of  these  will  now  be  obtained. 

Consider  the  inversion  formula.  As  ^  equation 

(3.17)  gives; 

<5o 

Z  TT  (bp'-fc")  as  Izl  and  Id  ~  (3.26) 

This  equation  states  that  Z  and  C  correspond  to  within 
a  real  scale  factor  at  large  distances  from  the  planing 
surface.  This  has  a  significant  meaning.  In  order  to 
obtain  the  outer  limit  of  the  inner  solution,  one  has  to 
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work  with  a  large  value  of 
value  of  ^  as  expressed  in 
Let 

A  -  ^0 

0  7,(bQ+CQ) 


Z  which  corresponds  to  the  large 
(3.26). 


(3.27) 


Far  away  frcan  the  planing  surface,  which  corresponds  to  the 
outer  limit  of  the  inner  problem,  equation  (3.17)  gives: 
(see  Appendix  A) 


Z  AQe'^^^f  (-CQ+t/l-CQ^)  5  +  logi;  (l+CQbjj-ibQ/l-CQ^) 

-  i  A-c +  ^^^log(bQ~^^^))  -  Cq 

-  (l+CQbQ)log(l+bQ)}  +  0(i)  (3.28a) 


Recall  that  log  C  is  treated  as  0(1)  as  approaches  a 
large  value.  To  the  first  term,  this  expression  can  be 
inverted,  as  follows: 


;  = 


A.(-c.+i/l-c.^) 


Z  +  0  (log  Z) 


(3.28b) 


Now  consider  the  one- term  inner  expansion.  Equation 
(3.8)  gives: 


‘o  = 


C+b, 


-]} 
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Far  away  from  the  planing  surface,  if  (3.28)  is  used,  this 
becomes : 

(-CqCos  a  +  A-CQ^sin  a)X  +  (c^sin  a  +  A-Cp^cos  o)Y 
+  0{log  X)  for  Icurge  ]z] 


In  the  outer  variables,  it  gives 


*^0^  {  “  CpCos 


+  (c.sin  a  +  /l-c-^cos 

^  n/i  x+hcot  a. 

+  O(log - - - ) 


where 


n  =  0{t)  =  0(e) 


The  one-term  outer  limit  of  the  one-term  inner  solution  is 

1  / - 3 

<5>0  ^  I(  -  CpCos  a  +/1-Cq  sin  a)x  +  (CpSin  a 

+ A-Cp^cos  a)yj  (3.29) 

Next  consider  the  inner  limit  of  the  outer  solution. 

The  one-term  outer  solution  is 


4>p  (x,  y,  z)  -v  X 


In  terms  of  inner  variables,  this  becomes: 
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(|>Q  eX  -  eHcot  a 

This  is  the  one- term  inner  expansion  of  the  one-term  outer 
solution.  To  be  matched  to  the  expression  in  (3.29) ,  it 
must  be  expressed  in  the  same  variables.  Thus  the  one-term 
inner  expansion  of  the  one-term  outer  expansion  is: 

<I>Q  (x,  y,  z)  'V.  X  (3.30) 

Equations  (3.25),  (3.29)  and  (3.30)  give: 

Cq  =  -  cos  a  (3.31) 

This  was  Icnown  to  Green;  however  he  assumed  it  to  be  true. 


IV.  SECOND-ORDER  SOLUTION 


IV.  1.  Second-Order  Outer  Solution 

The  boundary- value  problem  for  the  second-order  outer 
solution  can  be  obtained  from  equations  (2.13)  -  (2.16) 


(L)  =  0  in  fluid  (4.1) 


(FI)  •J'Oy'^ly  +  '*’oz'^lz^  ^"’l^^Oxy'^'Ox 


+  •J’a  4>a  +  ❖«  )  =  0  y=o 

Oyy^Oy  ^Ozy^Oz  ^ 


(4.2) 


”■2’  '*0x''lx  *  ^Ix^Ox  •  ♦ly  *  ♦02''lz  ♦lz''0x 


^  '^I'^Oxy^Ox  ”  ^l^Oyy  *  '^I'^Ozy^Oz  ® 

on  y=0 


(4.3) 


(R)  fix  =  ®  '  hy  =  hz 


0  X-^-w 

y-v— CD 


(4.4) 


If  the  equation  (3.24)  is  used,  the  free-surface 
boundary  conditions  become 


(FI')  (x,y,z)  =  0 


on  y=0 


(4.2’) 


^  ‘^Ox'^lx  *'  '^Ix'^Ox  ^ly  *  '^Iz'^Oz 

on  y=0 


(4.3') 
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The  equations  (4.1),  (4.2*) r  (4.3*)  and  (4.4)  constitute 
the  second  order  boundary  value  problem. 

The  equation  (4.2*)  states  that  the  velocity  potential 
can  be  continued  analytically  into  the  upper  half  space  as 
an  odd  function  with  respect  to  y=0.  The  equation  (4.1) 
states  that  tlie  potential  is  a  harmonic  function.  Up¬ 
stream  of  the  planing  surface,  the  potential  4)^^  can  be 
required  to  satisfy: 

(x,  0,  2)  =  0  on  y=0,  x<0,  (4.5) 

In  this  domain,  (|)^  must  be  continuous  at  y^O.  The  function 
is  thus  defined  in  the  entire  space  and  has  the  following 
property 


(|>j^  (x,  y,  z)  =  -  (x,  -y,  z)  (4.6) 

It  is  this  result  which  permits  the  outer  problem  of 
the  planing  surface  to  be  reduced  to  an  equivalent  airfoil 
problem.  Downstream  of  the  planing  surface,  equation  (4.2') 
then  becomes  just  the  usual  downstream  condition  for 
continuity  of  pressure  across  a  vortex  wake.  The  body 
condition  is  replaced  by  requiring  the  outer  expansion  to 
match  with  the  inner  problem.  The  line  of  singularities 
is  now  located  on  the  undisturbed  free  surface,  and  the 
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vortex  sheet  is  assumed  to  lie  on  the  x2-plane  which 
corresponds  to  the  free  surface  at  upstream  infinity. 

The  planing  problem  will  thus  be  solved  as  if  there 

were  a  vortex  wake.  Actually,  there  is  none;  instead  there 

is  a  free- surface.  It  is  simply  that  the  outer  problem  of 

the  planing  surface  is  mathematically  equivalent  to  a 

well-studied  airfoil  problem.  Concepts  such  as  circulation 

25 

have  no  place  in  the  planing  problem.  The  fluid  region 
is  simply  connected  and  the  potential  function  is  single 
valued.  There  is  no  circulation.  Nevertheless,  some  of 
the  terms  and  symbols  of  aerodynamics  will  be  used  in 
solving  the  problem. 

The  planing  surface  can  be  represented  by  a  line  of 
singularities  at  x=0  along  the  z-axis  in  the  unbounded 
fluid  as  e  approaches  zero.  The  strength  of  these 
singularities  is  still  unknown.  Since  there  exists  a 
lift  on  the  body,  one  must  provide  for  the  existence  of 
singularities  in  the  flow  behind  the  body.  These  are 
the  usual  wake  vortices,  which  can  be  represented  in 
terms  of  a  dipole  distribution  of  density  y{x,z)  in  tha 
xz-plane,  with  doublet  axes  normal  to  y=0  plane. 

The  wake  now  is  a  surface  of  discontinuity;  therefore 
the  distribution  of  doublets  must  be  extended  to  infinity 
over  the  wake  surface  W  as  well  as  over  the  planing 
surface  S.  The  span  of  the  planing  surface  is  located 
in  the  segment  between  z=-l,+l.  The  velocity  potential 
has  the  following  form 


_Lf‘- 

4ti  /  , 

J^llx 


2  2  , 
+y  +(z 


x>.  (c)d^ 

2^yS^U.„2,3/2 

(4.7) 


It  can  be  shown  that  the  equation  (4.7)  satisfies  the 
Laplace  equation  and  all  the  boundary  conditions  (4.1)  - 
(4.4).  The  first  integral  represents  a  plane  sheet  of 
dipoles  extending  to  infinity  downstream.  As  shown  in 
Appendix  B,  the  dipole  density  is  a  function  of  the  z 
variable  only,  y=y(z)-  The  second  and  third  integrals 
represent  discrete  lines  of  vertical  and  horizontal  dipoles, 
respectively.  Higher  order  singularities,  such  as 
quadripoles,  etc.,  need  not  be  considered  here  for  the 
second-order  approximation.  Sources  have  been  omitted,  on 
the  assumption  that  there  is  no  generation  of  fluid  in  the 
body. 


The  first  integral  can  be  integrated  with  respect  to  C 
to  give 


())j^(x,y,z) 


=  -i-T _ li 

+  JL^f  _ 

4n  r  2^  2^, 

+  A(C)di; 


[1  I  _ ^ _ _ 


Idr. 


,)2,3/2 


(4.8) 


44 


The  functions  y  (z) ,  X  (z) ,  and  p (z)  can  be  expanded  into 
series  which  are  asymptotically  valid  as  e  -►  0 


N 

y(z)  ^  I  Y„(z;t;) 
n=l  " 


> 


as  e  -*•  0 


u  (z) 


a. 


N 

I 

n=l 


lij^(z;£) 


as  e  0 


) 


X  (z) 


a. 


N 

I 

n=l 


X^(z;e) 


as  e  -*•  0 


(4.9) 


In  the  outer  problem,  as  aspect-ratio  AR  approaches 
infinity,  a  uniform  stream  is  recovered,  and  the  dipole 
density  eventually  vanishes.  Therefore,  one  has 


y^{z)  =  o  (1) 

=  o  (1) 

Xj^(z)  =  o  (1)  (4.10) 


Therefore,  the  two-term  outer  solution  is 


Y3^(C) 

~~T“  2 

y  +{z-0 


[1 


+ 


_ X 

[x^+y^+  (z-^) 


({)q  +  ({)q  X  + 
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+ 


+ 


[x^+y^+(z-C) 

tx^+y^+(z-c) 


2^72 


2j3/2 


dC 


dC 


(4.11) 


IV.  2.  Intezmediate  Step  Matching 

Let  !a=  1  and  n=2  in  the  equation  (2.6)  to  give: 

2-t  outer  expansion  of  (1-t  inner  expansion) 

=  1-t  inner  expansion  of  {2-t  outer  expansion) 

(4.12) 

In  order  to  find  the  former,  precede  as  follows: 

For  large  c,  the  inversion  formula  (3.17)  gives  a  two-term 
expansion: 


4  Of 

Z  'V,  AqC  +  AQ(e  -  bg)  log? 

+  AqO  -  tRjj)  +  0(-|-)  for  large  C  (4.13a) 

where 

Mq  =  Co  -  (l+boCo)  log(l+bo) 

Ro  =  bo/l-Co^  log  2  +  /l-Co^  log  (bo-*€^^~^) 

(4.13b) 


Next,  express  C  in  terms  of  Z  to  give 


46 


■  “o*  i  •  ‘"o  -  ^  °<5> 

U  «/ 

(4.13c) 

one-term  inner  solution  is  given  by  equation  (3.8).  Far 
away  from  the  planing  surface,  equation  (3-. 8)  gives: 

”0  "  -  ‘>o>  1'=%  - 

-  -0  ■  <‘>0'^0>  1°9<*>0'"=0> 

+  0(i)l} 

where  Aq,  Mq,  and  Rq  are  all  real  constants.  Next  take  the 
real  part  for  large  value  of  X  to  give 

<1>Q  -v  X  -  AgSina-tan”^  ^  “  ^0  ^ 

The  two-term  outer  limit  of  the  one-term  inner  expansion 
is  thus 

-1  Y 

-v  X  -  ApSina^tan  x  “  ^0  (4.14a) 

where 

Gq  =  Aq[MqCos  a  -  RQsina  -  cosa 


- (bQ-cosa)  log (bp-cosa) ] 


(4.14b) 
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Take  outer  variables 


^  xKhcQtoi  . 
u  e  0 


y-h  _ 
x+hcota  '’O 


X  tf  “1  V 

'»'  — {  X  ~  eA^sina  tan  x  ^  cota  -  cGq) 


A  constant  Kq ,  which  is  independent  of  space  coordinates , 
added  or  subtracted  from  a  velocity  potential  would  not 
violate  the  conditions  prescirbed  for  (Laplace 
equation  and  boundary  conditions) .  3o  the  two-term  cuter 
expansion  of  the  one-term  inner  expansion  is : 


i  (x  -  cA^sina  ♦  +  h  cota  -  eG^  -  eKjj) 


(4.15) 


The  two-term  outer  expansion  is  given  in  (4.11). 

To  take  the  irjner  limit,  this  solution  must  be  evaluated 
near  the  line  of  singularities.  Then  the  integrals  in 
(4.11)  become  divergent  integrals.  This  difficulty  is 
resolved  by  integration  by  part  to  give  (see  Appendix  C) 


(})q  +  x 


1 

(tan  tan 

z-C  X 


+ 


(z-C) 


7 


4tt 


x^+y^J_^  [x 


(z-C)  ^ 


dc 
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+  J _ ?L_  f  ! _ (2-0  .*^^1 


(4.16) 


The  strength  of  singularities  is  zero  outside  the  planing 
body,  which  gives: 


Y^(±l)  =  0  ;  Mjl(±1)  =  0  ;  Xj^(±l)  =  0 


Take  inner  varicibles 


(eX  -eH  cot  a)  -  ^(eY  +£H) 

4TT 


Yi(2) 


tan”^  (- 


eY+eH 


/ 

7  d^  2-! 
-'-I 


r) 


‘EX-eHcot  a' 

,  (eY+eH)m,(2)  +  (EX-EHcot  a)X,(Z) 

+  y-  '  . "■■— — ^  ^  —  I  .1  —  +  O(EY-i) 

(EX-EHcot  a)^  +  (eY+eH) 


(4.17) 


The  Cauchy  Principal  value  is  to  be  taken  for  the  integral. 

The  order  of  magnitude  of  Yj^tz)  with  respect  to  e  is 
not  known  yet.  However,  one  does  know  that  Y]^(z)  =  o(l). 
Therefore,  the  one-term  inner  expansion  of  the  two-term 
outer  expansion  is: 


(EX-cHcot  a)  - 


Yl  (Z) 

^  tan  (- 


eY+eH 


2it  'EX-EHcot  a 

^  (£Y+£H)Mj^(Z)  +  (EY-EHcOt  a)X^(Z) 

(EX-EHcot  a)^  +  (EY+EH) ^ 


) 
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In  the  outer  variables,  this  is  : 


From  the  equations  (4.12),  (4.15)  and  (4.18),  matching  gives 


Yj^(z)  =  2TreAQsin  a 

^^(z)  =  0 

Xj^tz)  =  0 

Gq  +  Kq  “  h  cot  a  =  0  (4.19) 

The  discrete  lines  of  vertical  and  horizontal  dipoles  w^(z) 
and  Xj^(z)  are  of  higher-order  than  that  of  So 

Y^(z)  =  0(e) 

M^(z)  =  Xj^(z)  =  o(e) 

4i^(x,y,z)  =  0(e)  (4.20) 


IV.  3.  Second-Order  Inner  Solution 

From  the  first-order  outer  solution,  5>q=0{~)  as 

expected.  However,  the  definition  of  the  asymptotic 

expansion  only  implies  that  $^=o(4'q).  The  exact  order  of 

2 

magnitude  of  is  still  unknown.  Therefore,  e  ^OZZ 
must  be  included  in  the  Laplace  equation. 


50 


Next,  use  the  s<me  argment  given  in  section  II.  4. 
the  possible  second-order  prc^lens  cd>tained  froai  the 
equations  (2.23)  -  (2.26)  are  included  :ln  the  folloving: 


(L) 

*lxx  *1YY  *  '  ® 

(4.21a) 

(FI) 

^♦ox*lX  *  "*0Y*1Y  *  “ 

(4.21b) 

(F2) 

♦ox«ix  *  ‘ix^ox  -  ^’lY  *  '  0 

(4.21c) 

(B) 

♦j^j^sin  a  +  ♦j^^cos  o  *  0 

(4.21d) 

Now,  the  order  of  will  be  estimated.  There  are 
two  possible  ways  to  do  this.  One  is  to  assign  every 
possible  order  of  magnitude  to  4^,  then  try  to  show  that 
each  assigned  order  of  magnitude,  except  one,  gives  rise 
to  a  contradition  among  the  Laplace  equation  and  all  the 
available  boundary  conditions.  The  exceptional  case 
which  is  acceptable  for  each  equation  is  thus  the 
possible  order  of  magnitude  for  The  second  approach 

is  to  consider  all  the  non-homogeneous  equations.  Determine 
the  order  of  magnitude  of  non-homogeneous  parts.  The 
lowest  one  will  be  the  possible  order  of  The  nice 

thing  about  the  method  of  matched  expansions  is  that,  if 
the  order  of  magnitude  is  not  correctly  estimated,  it 
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will  show  up  autowatically  in  the  later  part  durin9 
the  matching  performance.  %e  second  method  will  be 
used  here  to  estimate 

The  available  non-honcgeneous  equations  are 


0(e) 

*1XX  *1YY  -  ~  ^  ♦OZZ 


*1 

O(^) 


*1 

0(-^) 


0(1) 


^*ox*ix  ^*0Y^iy 


oz 


(4.22a) 


(4.22b) 


*0X**1X 


A  N  ^  £4  * 

*ix"ox  ^*iy 


(4.22c) 


Equation  (4.17)  gives: 

(4.22d) 

(4.22e) 

The  first  term  in  the  right  hand  side  of  (4.22e)  is 
matched  to  iiuier  solution  4q  and  the  second  term  is 
matched  to  iy  Therefore,  the  equation  (4.22e)  gives 


1- T.  I.  E.  of  2-T.  0.  E.  'V  0(ex) 

and 

2- T.  I.  E.  of  2-T.  0.  E.  'V  0(eX)  +  O(e^Y) 


'V 


0(e<I>Q) 


(4.22f) 
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The  kinematic  free-surface  condition  (4.22c)  vill  be 
used  to  determine  the  free-surface  deflection,  therefore,  it 
can  only  be  used  to  determine  the  order  of  From 

(4.22a),  (4.22b),  and  (4.22f),  the  lowest  order  for  4^  is 

4^  =  0(1)  (4.23) 

If  equation  (4.23)  is  used,  the  second-order  inner 
problem  becomess 

^IXX  ^  ^lYY  “  °  fluid  (4.24a) 

(FI)  2(4qjj4^jj  +  4^^$^^)  =0  on  Sj  (4.24b) 

^0X**1X  “^IX^^OX  “  ^*^1Y  “  ®  ®f  (4.24c) 

(B)  4j^jjSin  a  +  4j^y^os  a  =  0  on  body  (4.24d) 

Instead  of  this  boundary-value  problem  being  solved 
directly,  it  will  be  combined  with  the  first-order  inner 
problem.  The  new  boundary-value  problem  will  be  shown 
to  have  exactly  the  same  forms  as  the  first-order  inner 
problem,  and  the  solution  obtained  in  the  previous 


=  i 


i  ! 


f  5 

I,  ! 
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chapter  can  thus  be  carried  over  directly  to  give  the 
second' order  inner  solution. 

Consider  the  new  boundary* value  problem  in  which 
first-order  ana  second-order  conditions  are  added: 


*oxx  ’oyy  ixx  ’iyy 


in 


fluid  <4.25a) 


*0X  *0Y  ^^*0X*1X  *0y*lY^  “  ^  on  (4.25b) 


*0X**0X  “  ^*0Y  *  ^0X**1X  *  ^1X**0X  "  ^*1Y  "  °  (4.25c) 

on 


♦QjjSin  a  +  a  +  tj^j^sin  a  +  4>j^yCOS  a  =  0  (4.25d) 


on  body 


2  2 

Now,  consider  the  equation  (4.25b).  If  and 


are  added  to  the  left  side  of  equation  (4.25b),  the  effect 
will  not  be  felt  until  the  third  approximation  or  higher. 
Therefore,  without  loss  of  accuracy  to  the  second-order 


2  2 

approximation,  and  can  be  added  onto  (4.25b);  the 


same  argument  is  applied  on  (4.25b}  and  (4.25c).  One  obtains: 


(L) 

(fi  4.  A 

XX  YY 

(FI) 

*2  *2 

$  +  $ 

Y 

=  0 


in  fluid 


(4.26a) 


=  1 


on  S, 


(4.26b) 
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(F2)  **x  "  ^ 


on  S.  (4.2fc) 


(a) 


*  * 

4^sin  a  >  4^008  a 


=  0  on  body  (4.26d) 


where 


^0^  *1 


N  =  Nq  + 


(4.27) 


This  is  the  new  boundajry-value  problem.  The  governing 
differential  equation  is  the  two  dimensional  Laplace  equation. 
All  the  boundary  conditions  are  two  dimensional.  Therefore, 
the  solution  represents  a  plane  irrotational  flow.  The  method 

developed  in  the  first-order  inner  solution  can  be  applied 

* 

directly  to  this  new  problem.  Once  4  is  determined,  the 
second- order  potential  4^^  can  be  obtained  easily. 

In  the  first-order  problem,  a^,  b^  and  Cq  represent  the 
jet  thickness,  jet  direction  and  the  location  of  stagnation 
point  on  the  C-plane  respectively.  Now  due  to  the  three 
dimensional  effect,  these  values  will  be  changed  by  a  small 
amount . 

Let  a,  b,  and  c  represent  true  values,  assume  the  exist¬ 
ence  of  asymptotic  expansions  for  a,  b,  and  c 

N 

®  \io  ^n+r'  >®n> 

as  £-+0 


(4.28a) 
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"nio 


as  e-»0 


(4.28b) 


c 


N 

'^n^O  % 


c  ^-=o(c  ) 
n+1  n 

as  e-'-O 


(4.28c) 


For  the  new  boundary~value  problem  4^  ,  let 


* 

a 

* 

b 

* 

c 


ao  +  a^ 


‘>0  +  “l 


=0  *  'l 


(4.29) 


*  * 

Since  4  is  a  harmonic  fvinction,  ‘ts  conjugate,  ¥  ,  exists. 

* 

Let  the  analytic  function  F  be  the  complex  potential 


F  (z)  =  Fq(z)  +  Fj^(z)  =  4  (x,y)  +  i'F  (x,y) 


(4.30) 


Use  the  same  techniques  as  given  in  section  III.  1. 

* 

Map  the  F  -plane  onto  the  lower  half  of  ^-plane  as  shown  in 
figures  3.1  and  3.2  to  give; 

*  * 

F*(Z)  =  - F  [C  -  c*  -  (b*+c*)  log  -4^1  (4.31) 

IT  (b  +c  )  b  +c 

* 

Next,  define  the  new  mapping  function  to  be  H  (?;).  Then  the 
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inversion  formula  is: 


w  (Z(c' 


(4.32) 


where  W  (Z)  is  the  complex  velocity. 

i 

Next,  define  the  Kirchhoff  function  Q 


Q*  =  Qq  +  ^(0n+®i>  »  “  log  W*  (4.33) 


Then  on  the  free  surface  it  gives! 


1  -  =  0 


on  Sj 


(4.34a) 


On  the  body: 


Gq  +  0^^  -a 


0Q  +  Gj^  =  ir  -  a 


on  CB 


on  AC 


(4.34b) 


Now,  map  these  boundaries  onto  the  real  axis,  n^O,  of  the 
C-plane.  The  conditions  become: 


1  .  n<51)  = 


(n=0,|c|>l)  (4.35a) 


®0  +  ®1  =^TT-a 


(ti=0,Cq+Cj^  i<l) 
(n=o,-i<^<CQ+Cj^) 

(4.35b) 


0 


Equation  (4.35)  gives; 


W  (Z(?))  =  e 


ta 


C  -  (Cq+Cj) 


1  -  (Cq+Cj^)C  +  iAl-(cQ+c^)nu^‘ 


Then  the  inversion  formula  gives: 

*  -ia  *  *  *  r+b* 

z  =  <-c  (?-l)  ^  (l+c  b  )  log 

Tr(b  +c  )  b  +1 

log  (c+/c^-l) 

-  log  ub%  -  /(b;^-i)  (cizii) 

C+b 

Equations  (4.31)  and  (4.37)  provide  a  complete  second 
inner  solution. 


IV.  4.  Second  Order  Matching 

Let  ra=2  and  n=2  in  (2.6) ,  giving; 

2-T.  0.  E.  of  (2-T.  I.  E.) 

=  2-T.  I.  E.  of  (2-T.  0.  E.  ) 

For  large  C,  (4.37)  gives  two-term  expression. 


Z  =  c  +  A„  (e-<“ 


bp)  log  C 


t  AQe-^“(Mo  -  iRg) 


_  _VVfl^ 

TT  (bQ+C^  TT(b.+C„)^ 


1) 

(4.36) 


(4.37) 


■order 


(4.38) 


(4.39a) 


where 


(4.39b) 
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^nd  Aq,  Mq  ,  and  are  defined  in  (41. 3b). 
Next,  express  r,  in  terms  of  Z  to  two  terms: 


Z  _  Z  H  . 

—  —  —  I  +  t  ■  TT-  1 

Aq  Aq  Ajj  sin  a 


,  -id  w  \  1  Z 

(e  -  bj)  log  5- 


+  -  Mj) 


(4.39c) 


Substitute  (4.13b),  (4.29),  and  (4.39b)  into  (4.3i)  to  give: 


*0  ^  *1"  0^  t 

-  (e-*“-b(,)  log-l-  - 

-  (Cj+Cj)  -  (bj,+b^+C(,«^)  log-|- 

+  (bQ+Cg+bj^+Cj^)  log  (bg+bj^+Cg+Cj^) }  +  0(|) 

(4.40) 


Expressed  in  terms  of  outer  variables,  then  fhe  two-term  outer 
expansion  of  the  two-term  inner  expansion  is 


^0  ‘^l'^  e  sITTlt^  "  eAgsin  atan"-*-  J  hoot  a  -eGg-  eKg] 

(4.41) 


Now  consider  the  inner  limit  of  the  outer  solution.  The 
two-term  outer  solution  is  given  in  equation  (4.16).  It  is 
assiamed  that  the  derivative  of  Y]^(z)  with  respect  to  the 
space  coordinate  z  is  the  same  order  of  magnitude  as  y^(z). 
Then  the  two-term  inner  expansion  of  the  two-term  outer 


expansion  is: 
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-1  V  V 

^ '  -IF”  -i  -  h 

2^  ytij^(z)  +  xXj^(z) 


.r^dc. 

rj  dc  z-i 

-'-I 


2v 


2  2 
X  +Y 


(4.42) 


From  (4.38>,  (4.41)  and  (4.42),  matching  gives 


sin  g  f  d? 

4ir  J  dc  z-i 


(4.43) 


yj^(z)  =  Xj^(z)  =  0 


(4.44) 


The  equation  (4.44)  states  that  y^(z)  emd  X^(z)  are  of 
higher  order.  Therefore; 

Uj^(z)  =  Xj^(z)  =  o(e)  (4.45) 

From  (4.43),  depends  on  the  rate  of  change  of  Y2^(z) 
with  respect  to  the  spanwise  variable.  To  the  first  order, 
equation  (3.31)  gives 
a  =  cos  ^(“Cq) 

Now  the  hydrodynamic  angle  of  attack  will  be 

ttjj  =  cos“^(-Cq  -  c^) 

Therefore,  c^  represents  the  deviation  of  from  a. 

So  Cj^  represents  a  "downwash",  as  it  is  called  in  aerodynamics. 


V.  FREE- SURFACE  MATCHING 


V.  I.  Introduction 

The  planing-surface  problem  is  distinguished  from  the 
airfoil  problem  by  the  presence  of  a  free-surface.  In  Green's 
classic  two  dimensional  problem  of  a  plate  gliding  on  the 
surface  of  a  gravity-free,  infinitely  deep  fluid,  the  non¬ 
linear  free-surface  problem  can  be  solved  exactly.  However, 
Green's  solution  exhibits  an  anomaly.  The  free-surface 
elevation  is  given  by 

y  X  ~  log  |x|  as  |x|  *  (5.1) 

Thus  far  away  from  the  gliding  plate,  the  free-surface  drops 
off  logarithmically  to  ”«>.  The  location  of  the  gliding  plate 
can  not  be  prescribed  and  his  solution  is  not  unique. 

The  present  chapter  will  show  that  the  anomaly  exhibited 
in  this  classic  solution  is  just  part  of  a  proper  near-field 
expansion  of  the  complete  solution  which  includes  three 
dimensional  effects.  An  inner  expansion  does  not  necessarily 
satisfy  the  obvious  conditions  at  infinity;  it  must  only 
match  some  outer  expansion  in  a  proper  way.  A  far-field  free- 
surface  expansion  will  be  derived  and  matched  with  this  inner 
expansion.  Through  this  free  surface  matching,  the  height  of 
the  planing  surface  is  properly  defined.  The  depth  of  s\ib- 
mergence  can  thus  be  prescribed  and  the  solution  becomes  unique. 
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V.  2.  First-Order  Free- Surface  Matching 
A.  Near-Field  ^ree-Surface  Solution 

The  free  surface  and  planing  plate  on  the  physical  plane 
are  mapped  onto  the  real  axis  of  the  auxiliary  complex  c-plane. 
The  inversion  formula  (3.17)  gives  the  relationship  between 
the  physical  Z-plane  and  the  ?-plane.  To  describe  the  free^ 
surface  elevation,  let  n^O  and  on  the  equation  (3.17).  The 
image  of  the  free  surface  (n=0,|C|>l)  in  the  Cplane  can  now 
be  transformed  back  to  the  Z-plane: 


(5.2) 


X(C)  and  Y(C)  are  the  coordinates  of  a  point  in  the  free  surface. 
Far  away  from  the  planing  surface,  the  equation  (5.2)  gives: 


Z  =  X  +  iY  -  Aq(  C  +  (e’^^-bg)  log  C  +  e'^“(MQ-  iRQ)l+  0(i) 

(5.3) 

Aq,  Mq,  and  Rq  are  defined  in  (4.13b).  The  leading  edge  is 
mapped  onto  ?=-l  and  the  jet  onto  C'-bQ.  Therefore,  the  up¬ 
stream  free  surface  corresponds  to 


?  <  -bg  <  -1 

(5.4) 

log  ?  =  log  1  cl  -  -ni 

(5.5) 

Separate  the  real  and  imaginary  parts  in  (5.3)  and  eliminate 


the  parameter  ^  to  give  : 


Y  “v  -  A^sin  a  logjxi  +  Anisin  a  log  Aq  -  7r(cos  a  -b^) 

-  M^sin  a  -  RqCos  a]  +  O(log-L^)  (5.6) 

This  is  the  free-surface  elevation  far  upstream  and  dovmstream 
frcHn  the  planing  surface.  As  X  ia>,  the  free-surface  deflection 
beccmes  Y  ■*  -<*>.  This  anonaly  is  also  observed  in  Wagner's  and 
Green's  two  dimensional  solutions.  However,  it  will  be  shown 
in  the  present  work  that  the  equation  (5.6)  is  just  a  proper 
inner  description. 

B.  Par-Field  Free-Surface  Solution 
Equation  (3.20)  gives 

^  ♦oz''0z  =  "  °”  y'® 

If  (3.24)  is  used,  equation  (5.8)  becomes 

"ox  "  ® 

Integration  gives 

Hq  =  0(2)  (5.8) 

where  D  is  function  of  z  only.  However  (3.24)  represents  a 
uniform  stream  and,  to  the  same  order  of  magnitude,  the  free- 
surface  will  correspond  to  the  undisturbed  free-surface  elevation. 
This  gives  D=0.  So 


0 


(5.9) 
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This  will  be  called  the  "one-term  expansion"  of  the  free-surface 
elevation . 

C.  First-Order  Free-Surface  Matching 
Equation  (3.29)  gives 

(  -CqCos  a  +  A-CQ^sin  o)x  +  (CQsin  o+A-Cq^cos  a)y 
X  as  Izl  "o  (5.10) 

where  c^  =  -  cos  a  is  used.  This  is  the  one-term  outer  expansion 
of  the  one-term  inner  expansion  of  the  velocity  potential,  which 
corresponds  to  a  uniform  stream  at  its  outer  limit.  Feu:  away 
the  free-surface  can  thus  be  represented  by 

Nq  =  0  on  Sj  (5.11) 

This  is  the  one-term  outer  expansion  of  the  one-term  inner 
expansion  of  the  free-surface  elevation. 

From  (5.9),  ta)ce  inner  variables  to  obtain  one-terra 
expression  and  then  converted  back  to  outer  variables  to  give 

tIq  =  0  (5.12) 

This  is  the  one-term  inner  expansion  of  the  one-term  outer 
expansion. 

Equations  (5.12)  and  (5.13)  agree  with  the  matching 
principle  (2.6).  The  first-order  free-surface  matching  does 
not  provide  additional  information;  in  a  sense,  it  is  trivial. 
Nevertheless  this  matching  does  establish  the  relationship 
so  that  the  higher-order  free-surface  expansions  will  follow 
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the  matching  principle. 


V.  3.  Second-Order  Outer  Free-Surface  Solution 

Equation  (4.3*)  gives  the  kinematic  free-surface  condition: 


♦ox''lx  *  ♦ix^Ox"  *1,  +  hz^Oz  =  0  on  y=0  <5-«> 


If  (3.24)  and  (5.9)  are  used,  one  obtains: 


Integration  with  respect  to  x  gives 


^"1  y»0 


0,  z)dC 


(5.14) 


(5.15) 


The  lower  limit  is  taken  to  be  -»  to  coincide  with  the 
radiation  condition.  In  appendix  D,  it  is  shown  that 


where  the  Cauchy  principal  value  is  taken.  Substitute  (5.16) 
into  (5.15)  and  change  the  order  of  integration  to  give 


The  two-term  outer  solution  is  thus 
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(5.18) 


V.  4.  Intermediate  Free-Surface  Matching 
Let  m=l  and  n=2  in  (2.6)  to  give 
1-T.  I.  E.  of  (2-T.  0.  E.) 

■  2-T.  0.  E.  of  (1-T.  I  .  E.)  (5.19) 

Siibstitute  inner  variables  in  (5.18)  to  give 

no+  “  |^Yi(*Ul+log(ex-eHcot  a|l 

2|Z-c|dC 

+  o(e^)  (5.20) 

C>0 

C<0  (5.21) 


where 


w.; 

'C 


sgn  (Z-c)  »  { 


dYi (C) 

— — sgn  (Z-C)  log 

^^1^^^  cX-cHcot  a., 

+1 
-1 


z- 

z- 


So  the  one -term  inner  expansion  of  the  two- term  outer 
expansion  is: 


^0  '’l'^  ■  log|x|] 


+ 


dYi (C) 
-dc - 


(2-C) 


log  2|z-i;|dc 


(5.22) 


Next,  substitute  outer  variables  in  (5.6)  to  give  the  two 
term  outer  expansion  of  the  one-term  inner  expansion: 
y  'v  h  -  eApSin  alog  jx]  +  eA^sin  alog  e+  eE^ 


(5.23) 


O 


F 
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where 

Eq  ~  AgTsin  alog  A^  -  ir(cos  oi-bQ)  -  N^sin  o-  R^cos  o] 

(5.24) 

f'rom  (5.19),  (5.22)  and  (5.23),  the  intermediate  matching  gives 


h  =  -  elog  eA^sin  a  -  eE^  -  Yj(z) 


i_r 

4TrJ  d? 


-sgn  (z-O  log  2|z-c|dc 


(5.25) 


where  is  gxven  in  (5.24) 

The  leading  term  in  the  inner  solution  is  -cA^sin  alog|x|; 
this  term  is  exactly  matched  to  one  term  in  the  outer  solution 
-  ^  Y;i^(z)  log  |x|.  This  shows  that  the  anomaly  in  the  classic 
solution  is  just  part  of  a~proper  near-field  expansion  of  the 
complete  solution.  It  is  also  noted  that  far  away  from  the 
planing  surface,  if  AR  approaches  infinity  or  e->-0,  physically 
it  corresponds  to  a  large  span  of  small  chord  length  and  the 
planing  surface  shrinks  to  a  line.  To  this  lowest  order, 
there  is  no  disturbance  in  the  far  field,  which  agrees  with 
the  description  in  (5.25). 

The  location  of  the  planing  surface  is  now  uniquely 
defined  with  respect  to  the  undisturbed  free  surface.  In 
natural  variables, (5.25)  gives 
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vrtiere  are  defined  in  (4.13b).  This  is  the  important 

result  of  this  chapter. 


V.  5.  Second-^Order  Inner  Free-Surface  Solution 
Far  away,  the  inversion  formula  gives: 

Z  =  (Aj  +  Aj  +  +  A„(e-<“-bo)  log  C 

+  Aje‘‘“(Mj-iRj)  +  0(i)  (5.27) 

On  the  free  surface: 


Z  =  X  +  tY 


(Ao+Ai+iAo5^)C  +  Ao(e-i“-b(,)log  5 

+  v‘'“<v«o>  ^  °'r’  as  Id- 


Separate  the  real  and  imaginary  parts  and  eliminate  the 
parameter  C  to  give; 


y  'V/ 


sin  a 


-  ApSin  alog 


X|  +  E.  +  0( 


logj  X.I 

X 


(5.29) 


V. 


6 .  Second 
Let  m=2 


Order  Free-Surface  Matching 
and  n=2  on  (2.6) ; 


2-T.  I.  E.  of  2-T.O  .  E. 

=  2-T.  0.  E.  of  2-T.  I.  E.  (5.30) 

From  (5.20),  the  two-term  inner  expansion  of  the  two-term 
outer  expansion  is 

Yi(2)  f  dr  (c) 

•^0  +  ‘  “27”  I"'!'  ^.47>i— 3c - 

,  ,  1  P  (C) 

•  2h-c|dC  -  ^  dC  (5.31) 

In  (5.29),  use  outer  varifibles  to  give  the  two-term  outer 
expansion  of  the  two-term  inner  expansion  : 

c. 

y  tiQ  +  h  +  -  eApSin  alog  |x| 

+  elog  e  A^sin  a  +  eA^Isin  alog  Aq 
-  Tr(cos  ot-bp)  -  MpSin  a  -  R^cos  a]  (5.32) 

From  (5.30),  (5.31)  and  (5.32),  matching  show  that  each  term 
in  (5.31)  is  exactly  matched  by  one  term  in  (5.32),  and  vice 
versa.  In  a  sense,  this  matching  is  trivial.  However,  it 
does  give  a  hint  that  so  far  the  analysis  is  correct. 


VI.  HYDRODYNAMIC-LIFT  CALCULATION 


VI.  1.  Chord  Length  of  the  Planing  Surface 

The  hydrodynamic  problem  can  be  completely  described  by 
parameters  a^,  b^,  and  Cq  for  the  first-order  solution  and  a^^, 
bj^,  and  c^^  for  the  second-order  solution  respectively.  However, 
thoje  parameters  are  related  to  the  chord  length  at  each  section. 
Let  the  chord  length  be  i(z)  in  natural  varicd>le  and  L(z)  in 
inner  variable.  The  chord  length  is  mapped  onto  the  C'pl^ne 
by  the  segment  |C|  <  1  on  the  real  axis  in  Figure  3.2. 

The  trailing  edge  and  leading  edge  of  the  plate  are  mapped 
onto  1»-1  respectively.  From  equation  (3.11) ,  the  chord 
length  is  found  to  be 


L(z)e‘^“  =  “  Z 


=r  H(c' 

I  W(Z(c' 
-'-I 


T) 


where  and  Z.  denote  the  trailing  edge  and  leading  edge 

i  ii 


locations ,  respectively .  So 


i(z)  =  Ee 


H(C' 

J  ^  WTTT 


L 


dC 


(6.1) 


A.  First-Order  Chord  Length  Expression 

If  (3.11)  and  (3.18)  are  substituted  into  (6.1),  it 
follows  that: 
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e6, 


bn-l 


+  lT»4^^(bQ-^^)l 


(6.2) 


The  first-order  parameters  a^,  b^,  and  Cq  are  thus  related  to 
the  chord  length  Hz). 

B.  Second-Order  Chord  Length  Expression 

Substitute  (4.32a)  and  (4.36)  into  (6.1)  to  give 


£ 


e(6-+6T)  b^+b,-l 

IT  (bQ+bj_+CQ+Cj^)  (l+(Co+Ci)  (bo'‘’^l\  bp+bj^+l 

-2(Cq+Cj^)  +  tt/i- (Cq+Cj^)  ^ (bQ+b3^-/(bQ+bj^) ^-1)  1  (6.3) 


Thus,  the  parameters  a^^,  b^^,  and  are  also  related  to  the 
chord  length  £(z). 

VI.  2.  Wetted  Length  of  the  Planing  Surface 

A  jet  is  thrown  off  at  the  leading  edge.  Therefore  it 
is  somewhat  ambiguous  to  refer  to  the  "wetted  length"  in  this 
case.  However,  usually  the  jet  thickness  is  small  compared 
to  the  chord  length.  The  "pressure  length  can  thus  be  defined 
as  the  distance  between  the  trailing  edge  and  the  stagnation 
point  on  the  planing  plate.  Then  define  the  "wetted  length" 
L^(z)  as  the  pressure  length,  i.e., 
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where  denotes  the  stagnation  point. 

To  the  first-order,  it  is 


‘w<-'  =  7(B^  ^  -  <=0  ^  <=o' 

+  +  boSln-l(-Cg) 


0  —1 
+  /b  -1  sin 

0  0 


(6.4) 


To  the  second  order,  it  is 


e(«Q+6j^) 


bo+bi+1 


'  ir(bQ+bj^+CQ+Cj^)  <l+(Co+Ci) 

-  (Cq+c^)  +  (Cq+Cj)^  +  /l-  (Cq+c^)  ^  I j(bo+b^-/ (bQ -b^)  ^-1) 


-  A-  (Cq+Cj^)  ^  +  (bQ+bj^)sin"^(“CQ-Cj^) 

+  /(b,,4b,)2-l.sin 


'0  1' 


(6.5) 


VI.  3.  Jet  Direction 

Let  jet  angle  3  be  the  angle  the  jet  makes  asymptotically 
with  the  plate.  The  jet  J  in  the  Z-plane  is  mapped  onto  ^=”bQ. 
Equation  (3.18)  gives 


,  ...  .  l+b^Crt-i/l-c-^ /hJ^-1 

(-Tr+a-3)^_  gta _ 0  0  0  0 


-b,^-c 


The  real  part  gives 
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„  -l/  0*^0, 

£  =  cos  (  --r . -r-—) 

^0  ^G 


(6.6) 


VI.  4.  Lift  Force 

Lift  force  can  be  obtained  by  integrating  pressure  along 
the  plate.  This  is  a  near-field  problem. 

Let  P  be  the  pressure  in  the  fluid,  F  be  the  force  and  t 
be  the  lift  on  the  planing  surface  per-unit— length  in  spanwise 
direction.  For  a  steady,  irrotational,  incompressible  flow, 
the  Bernoulli  equation  gives 


P 

P 


+ 


+  ) 


(6.7) 


where  p  is  the  fluid  density.  Let  P^=0  and  use  inner  variables 
to  obtain: 


P 

P 


+  -(♦  ^ 
2'  X 


2  2  2 
+  \  ) 


(6.8) 


Now  assume  the  following  asymptotic  expansions  exist  for  P,F, 
and 


N 

P  'V'  Z-  P 

n=0  n 


N 

F  -v  F 
n=0  n 


as  e-^0  for  fixed  X  ,y 

as  e-*-0  for  fixed  X 


73 


T 


T 


n 


<«•»> 

for  fixed  X,Y 


Substitute  (6.9)  into  (6.8)  to  give 


-r  *  T‘*0x^  ^  ‘OY^’  *  -r  ♦  ♦ox*lX  "■  *0Y*1Y 
^  r  *  l'*lx"  ^  ”0X*2X  ♦  ^lY^  ^  "0Y*1Y  ^ 


1 

I 


(6.10) 


A.  First-Order  Lift  Force 

Keeping  only  the  leading-order  terms  in  (6.10),  one  obtains 
^0  ^1,4  2  .  *  2.  1 

T  *  *0y  ^  2  (6.11) 

The  force  on  the  planing  plate  is 
rt.e 

Fq  =/  PgdZ  (6.12) 

J  l.B. 

where  Z.e.  and  t.e.  denote  the  leading  edge  and  trailing  edge 
of  the  planing  plate  respectively.  The  integral  is  taken  along 
the  planing  plate.  This  force  is  normal  to  the  plate  and  can 
be  decomposed  into  vertical  and  horizontal  components,  that  is, 
lift  and  drag,  respectively.  The  lift  force  per- unit-length 
can  now  be  obtained  from  equation  (6.12)  to  give: 
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rt.e. 

Tq(z)  =  =  ffg  1 1  |(l-lWjj|^)dZ)  (6.13) 

Jz.e. 


If  equations  (3.16)  and  (3.10)  are  substituted  into  equation 
(6.13),  one  obtains: 


(bQ-^^^)  *^^^cos  a 


(6.14) 


where  a^,  bQ,and  can  be  obtained  by  solving  three  simul¬ 
taneous  equations  (3.31),  (5.26)  and  (6.2).  The  equation  (6.14) 
gives  the  lift  force  per-> unit— length  in  spanwise  direction. 

This  equation  agrees  with  Green's  non-linear  solution.  However, 
Green  could  only  compute  the  lift  force  in  two  special  cases. 

One  case  is  bQ-*-l,  which  corresponds  to  the  infinitely  long 
planing  plate,  that  is,  the  very  lightly  loaded  planing  plate. 
The  other  case  is  bQ'^^r  corresponding  to  Rayleigh's  cavity 
flow  which  is  of  no  interest  in  the  planing  problem.  In  the 
present  work,  the  lift  force  can  be  computed  uniquely  for 
arbitrary  angle  of  incidence  and  for  a  heavily- loaded  as  well 
as  a  lightly- loaded  planing  arface.  Numerical  calculations 
will  be  presented  in  the  next  section. 

B.  Second- Order  Lift  Force 
Equation  (6.10)  gives 
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C  • 


5 


r 


t 

f 


The  force  per-unit— length  on  the  planing  plate  to  the  second- 
order  is: 


^0  *  ^1 


J-t  <e. 

(Po«l)dZ 

X*  •  6  • 


(6.16) 


If  equations  (4.32),  (4.36)  and  (6.15)  are  substituted  into 
equation  (6.16),  the  lift  force  to  the  second-order  is  found 


to  be 

Tq(2)  +  Tj^(z) 


®0^^1  >C 

(bo+bj-/(b^^ 


-(Cfl+Ci) 


-l)cos  a 


2 


(6.17) 


VI.  5,  First-Order  Numerical  Calculation 

There  are  three  unknowns  2md  three  available  equations 
(3.31),  (5.26),  and  (6.2).  Therefore,  it  appears  likely  that 
a  solution  can  be  found.  It  is  noted  that  c^  is  function  of 
a  only  ^md  is  easily  obtained.  The  problem  can  thus  be  reduced 
to  two  unknowns  with  two  equations  as  follow: 


h  = 


^0  ^0 


-  MgSin  a  -  Rq  cos  o+  sin  a]  + 


sgn  (z-^)*log  2lz-c|dc 


sin  af  fd  .  ^0 

2^  J_1  ^ 


-) 


0 

(6.18) 


=  n  (bjc^l  ^  ^  -  2Cq 

+  Ti/i-CQ^(bQ-^Q^-l)] 


(6.19) 


r 


E 
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Now  can  be  expressed  as  a  function  of  b^,  froa  (6.19).  If 
a^,  sc  obtained,  is  substituted  back  into  (6.18),  it  gives 
a  non-linear  integral  equation  of  one  variable  b^  with  a 
complicated  kernel.  To  perform  a  numerical  calculation,  it 
is  more  convenient  to  solve  a  system  of  two  non-linear 
equations  (6.18)  and  (6.19)  simultaneously  than  a  non-linear 
integral  equation  with  a  very  complicated  kernel. 

Now,  express  (6.18)  and  (6.19)  symbolically  as  following: 

h  *  g(aQ*”^  bQ^”h  +  I.P.(aQ^''“^^ 

l(z)  =  (6,20) 


Where  I.P.  represents  the  integral  term  in  (6.18),  and  the 
superscript  '  ^  denotes  a  value  at  the  n-th  iteration.  For 
the  first  iteration,  integral  part  is  set  equal  to  zero. 

The  problem  requires  the  solution  of  the  two  non-linear 
algebraic  equations.  Then  the  iteration  is  terminated  if 


^  (n)_  ^  (n-1) 


THT 


<  6 


and 


.  (n)  ,  (n-1) 

*’0  ‘’o 


(n) 


<  6 


(6.21) 


where  6  is  some  positive  number  which  is  much  less  than  one. 
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A.  Cusp  Shape 

Let  the  chord-length  distribution  be 

l(z)  =  (6.22) 

fdiere  is  the  chord  length  at  middle  section,  z=0. 

The  chord  length  £(z)  and  its  first  derivative  with 
respect  to  z  are  continuous  at  the  tip,  as  elsewhere.  Recall 
that  h  is  the  height  of  the  trailing  edge  above  the  un¬ 
disturbed  free  surface;  physically,  it  can  be  related  to  a 
loading  condition  on  the  pleming  surface.  If  the  loading 
is  heavy,  the  trailing  edge  will  be  located  much  below  the 
tuidisturbed  free  surface  and  consequently  h  will  have  a 
large  negative  value. 

For  £q=0,4,  0=15®,  the  coefficient  of  lift-per-unit  • 
span  (local  lift  coefficient)  formed  with  the  wetted  length 
is  given  in  Figure  6.1.  for  different  value  of  h.  Note 
that  the  ordinate  scale  does  not  start  at  zero.  The  jet 
thickness  Sq  at  each  section  is  shown  in  Figure  6.2. 

A  cusp  shape  is  not  a  practical  shape.  More  detail 
will  be  discussed  for  the  elliptic  shape. 

B.  Elliptic  Shape 

Let  the  chord-length  distribution  be 

i(z)  =  fpd-z^)®*^  (6.23) 

The  chord  length  is  continuous  at  the  tip,  as  elsewhere; 
however  its  first  derivative  is  not  continous  at  the  tip. 
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Observation  frm  experiments  indicates  that  tip  effects  are 
local  phenomena,  and  the  solution  obtained  from  the  present 
theory  can  be  expected  to  give  a  good  result  except  in  the 
very  neighborhood  of  the  tips. 

Local  lift  coefficient  distribution  in  spanwise  direc¬ 
tion  for  a=15®,  10^,  and  5®  are  plotted  in  Figures  6.3,  6.4, 
and  6.5  respectively  for  the  case  4^=0. 4.  Frtxn  the  defini¬ 
tion  of  b^,  the  value  of  bQ  is  bounded  by 


For  the  case  if  equations  (6.2)  and  (6.4)  are 

used,  one  obtains: 

where  4^  is  the  chord  length  in  natural  variable.  From 
equations  (6.4)  and  (6.14),  numerical  calculation  gives: 

0.289  for  0=15° 

— 0.214  a=lQ° 

pU‘^4„  0.120  a=5  (6.23) 

w 


The  upper-bound  limit  for  the  local-lift-coefficient  is 
thus  obtained.  Note  that  this  upper-bound  limit  is  a 
function  of  a  only.  This  statement  can  also  be  seen  in 
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in  Figure  6,8. 

For  the  case  bQ-*^,  if  equations  (6.2)  and  (6.4)  are 
used,  one  obtains: 

y-  -v  1  as  b^-^" 

Now,  if  equation  (6.6)  is  used,  the  jet  angle  8  is  found 
to  be: 


B  180®  as  b^-^ 

This  is  a  cavity-flow.  However,  this  does  provide  a  lower- 
bound  for  the  extremely  heavily  loaded  planing  surface. 
Numerical  calculation  gives: 

T.  0.167  for  a*15® 

-y”-  -v  0.121  a=ip° 

pU^Jl^  0.066  0=5°  (6.24) 

If  equation  (6.18)  is  used,  the  depth  of  submergence  h  is 
found  to  be 


h 


as  bp-^® 


Those  upper-bound  and  lower-bound  limits  are  plotted 
in  Figures  6.3,  6.4,  and  6.5  respectively.  Note  that  due 


f 

I 

f 
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Figure  6.9  Local  Lift  Coefficient  vs 
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to  the  three-dimensional  effect,  the  values  of  upper-bound 
and  lower-bound  local  lift  coefficients  decrease  along  the 
spanwise  direction  and  approach  zero  at  the  tips. 

The  local  lift  coefficient  at  the  middle  section  is 
plotted  with  respect  to  angle  of  attack  a  in  Figure  6.6. 

In  Figure  6.7,  the  ratio  of  wetted  length  to  chord 
length  »  at  the  middle  section  is  plotted  against  the 

height  of  the  trailing  edge  for  different  values  of  ^q. 

In  Figure  6.8,  at  the  middle  section,  is  plotted 
against  local  lift  coefficient  for  different  values  of 
for  a=15^.  It  shows  that  for  a  small  value  of  say, 

fjg<0.15,  the  chord  length  has  no  significance.  The  hydro- 
dynamic  results  are  the  same  whether  is  0.8,  0.6,  0.4, 
or  0.2. 

In  Figure  6.9,  is  plotted  against  the  local  lift 
coefficient.  This  figure  can  also  be  obtained  from  Green's 
non-linear  solution.  It  is  seen  that  the  local  lift  coeffi¬ 
cient  takes  a  constant  value  for  a  fixed  a  when  is  small. 
For  the  case  a-»‘0,  from  equations  (6.4)  and  (6.14),  one  obtains 

^0  1 

— ^  -v  Ua  as  a-^0  (6.25) 

pu\ 

Wagner  (1932)  considered  the  linearized  case  of  an 
infinitely  long  plate  and  very  small  a.  Wagner  gave  the 
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result  as  given  in  equation  (6.25).  Therefore,  the  present 
non-linear  theory  does  approach  Wagner's  linearized  solution 
as  of'D.  For  a  finite  value  of  o,  the  present  theory  gives 
the  values  of  local  lift  coefficient  in  equation  (6.23) , 
while  Wagner's  linearized  theory  gives  the  value  of  0.411, 
0.274,  and  0.137  for  a=15®,  10®,  and  5°  respectively.  It 
is  noted  that  a  non-linear  effect  roust  be  considered  if  a 
is  not  very  sroall. 

VI.  6.  Second-Order  Nuroerical  Calculation 

In  order  to  obtain  a  unique  second-order  solution,  the 
parameters  a^,  b^,  and  c^  roust  be  determined.  There  are  three 
vinknowns,  but  so  far  only  two  equations  (4.45)  and  (6.2b), 
have  been  obtained.  One  is  obtained  from  the  velocity  poten¬ 
tial  matching,  and  the  other  is  from  the  chord-length  rela¬ 
tionship.  The  second-order  free-surface  matching  fails  to 
provide  an  additional  equation  to  define  the  depth  of  sub¬ 
mergence  . 

From  the  first-order  solution,  it  appeared  that  a 
second-order  outer  solution  had  to  be  found  to  provide  a 
height  reference  for  the  first-order  complete  solution. 
Accordingly,  a  third-order  outer  solution  must  be  found  in 
order  to  provide  a  height  reference  for  the  second-order 
complete  solution.  An  equation  similar  to  (5.26)  will  then 
be  obtained  from  the  free  surface-matching.  Then  the 


/ 
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second-order  hydrodyn^Ulic  problon  will  be  completely  solved. 
The  third-order  outer  solution  will  be  formulated  in  the 
next  chapter.  However,  at  the  present  stage,  two  special 
cases  can  be  readily  calculated  before  entering  the  complicat¬ 
ed  third-order  outer  solution. 

A.  b  =  b  +  b-  1 

o  1 

This  corresponds  to  a  lightly  loaded  planing  surfaces. 

The  solution  similar  to  Wagner's,  with  a  second-order 
correction,  can  be  obtained. 

B.  Angle  of  Attack  a  Is  Small 

The  parameter  c^  vdiich  represents  a  downwash  correction, 
is  given  in  (4.43).  It  is  a  function  of  a^,  b^,  Cq,  and  a 
only  and  can  be  obtained  easily  from  the  first-order  solution. 
Now,  if  the  angle  of  incidence  is  small,  the  problem  can  be 
linearized  as  in  aerodynamics,  so  that  the  second-order  effect 
is  mainly  due  to  downwash.  One  c^tains: 

T„(z)  +  0  ^  COS  a 

°  °  ^  (6.26) 

The  solution  of  Cj^  is  given  in  Figure  6.10  for  the  case 
of  an  elliptical  planing  surface.  As  expected,  the  solution 
for  Cj^  breaks  down  at  the  tips.  However,  the  curves  indi¬ 
cate  that  it  is  quite  a  weak  singularity  at  the  tips.  The 
solution  can  be  expected  to  be  valid  all  across  the  span 
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except  in  the  very  neighborhood  of  the  tips.  Local  lift 
coefficient  is  given  in  Figure  6.11.  For  a  heavily  loaded 
planing  surface,  the  downwash  correction  is  more  significant 
than  a  lightly  loaded  planing  surface. 


VII.  THIRD-ORDER  OUTER  SOLUTION 


VII.  1.  Introduction 

In  order  to  obtain  a  complete  second-order  solution,  it 
is  necessary  to  solve  the  far-field  problem  to  three  terms. 

A  third-order  free-surface  expression  relating  parameters  a^^, 
bj^,  and  c^^  to  the  depth  of  submergence  similar  to  (5.26),  can 
chen  be  obtained  from  tlie  matching  of  the  three-term  outer 
expansion  of  the  two-term  inner  expansion  to  the  two-term 
inner  expansion  of  the  the  three-term  outer  expansion.  The 
second-order  problem  can  thus  be  c(xnpletely  solved. 

VII.  2.  Third-Order  Outer  Boundary- Value  Problem 

From  (2.13)  -  (2.16),  the  following  are  obtained: 
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‘*’2xx 

+ 

* 

^2zz 

=  0 

in 

fluid 

(7. 
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.4) 

The  free-surface  boundary  condition  (FI)  can  be  rewritten 
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♦2X  '  -  I-""  <’-5> 

The  right-hand  side  of  the  above  equation  is  known  from 
the  first-and  second-order  outer  solutions.  Therefore, 
combined  with  Laplace  equation  (7.1),  the  problem  is  a  non- 
hcMnogeneous  boundary-value  problem.  The  governing  differen¬ 
tial  equation  is  linear  and  so  the  principle  of  superposition 
is  applicable.  The  final  solution  can  be  obtained  by  adding 
the  homegeneous  solution  to  the  particular  solution. 

Define  the  velocity  potential  ®  particular 

* 

solution  of  this  non-homogeneous  problem.  And  let  p  be 


* 

F 


»  -  £  (A  2  4^ 

2  '^Ix  ^ly 


2+ 


on  y=0 


(7.6) 


Then  the  dynamic  condition  (Fl)  becomes 


2x 


(7.7) 


The  equation  (7.7)  can  be  interpreted  in  terms  of  a  pressure 

distribution  over  the  free  surface.  This  problem  will  be 

solved  by  the  use  of  a  double  Fourier  transform.  (See  Wehausen 
31 

&  Lai tone  ) 


VII.  3.  Third-Order  Outer  Solution 
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Define  the  double  Fourier  transfoirm  with  respect  to  x 
and  2  to  be 


dx*dz 


(7.8a) 


Then  the  inverse  transform  is 


(^(x,y,z) 


■^ri 


“'de-dk 

(7.8b) 


Now  take  tlie  double  Fourier  transform  of  the  dynamic  condition 
(FI)  to  give: 


sec  6  ;v* 
tkp  ^ 


on  y*0  (7.9) 


Next^take  the  transform  of  the  Laplace  equation  to  give 


0 


This  is  an  ordinary  differential  equation.  A  general 
solution  is 

^2p(Jt/y/9)  =  A(k,y,0)e~^^  +  B(k,y,e)e^^  (7.10) 

Note  that  y  vertically  up  is  positive.  If  radiation 
condition  is  used,  equation  (7,10)  becomes: 
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=  B{k,y,e)e^y 


(7.11) 


From  (7.9)  the  constant  B  it.  determined  to  give 


t>  #1.  ..  a\  _  sec  6 

■  ^kT  P 


So  the  equation  (7.11)  gives 


(7.12) 


^00  ^  00 

sec  0  /  I  * . 


^)^-fk(xcos  e+zsin 

(7.13) 


Take  the  inverse  transform  to  give 


♦ip** 


^Tt  —00  -oo 

"  hrjoj-z  y__^|^ycn,y,nl 

.gik(incos  e+nsin 


,gik(xcos  0+zsin  6)j3]^.j30 


Integration  with  respect  to  k  gives; 


<l>2p{x,y,z) 


4tj  J-v/2  J-«> 

(x-m)cos  0+(z-n)sin  0 

*  2  2 
y  +[(x-m)cos  0+(z-n)sin  0] 


d0 ’dm’dn 

(7.14) 


where  (tij^(x,y,z)  and  ^^^(XyZ)  are  given  in  (4.11)  and  (5.17). 
The  strength  (z)  (z)  »  and  Xj^(z)  are  given  in  (4.19). 


94 


Equation  (7.14)  gives  a  particular  solution.  The 
homogeneous  solution  is  obtained  frcwi  (4.8);  it  gives: 


1  yf,*?) 

•y ~  471/  2  2^^ 

J-1  y  +(z-o^ 

1  fi 

4Tr  .  TTm^ 
J-1  [x  +y  +(z- 

1 

^V-1  [x2+yV(z- 


^)^  /x‘^+y‘‘+(z-?) 

yy2(^) 

2  2 
+  (z-0  1  ' 


?)V7I 


(7.15) 


The  first  integral  in  the  above  equation  represents  the 
correction  >'2^2)  to  the  vortex  strength  Yj(z)  in  (4.9). 
The  final  third-order  outer  solution  is 


<t»2(x,y,z)  =  c|)2p(x,y,z)  +  <{>2h(x,y,z)  (7.16) 

The  unknowns  ^2^^^  ^  determined  by 

intermediate  matching  to  the  second-order  inner  solution. 


VII.  4.  Intermediate  Velocity-Potential  Matching 
Let  m=2  and  n--3  in  the  equation  (2.6); 

2-T.  I.  E.  of  3-T.  O.  E. 

=  3-T.  0.  E.  of  2-T.  I.  E.  (7.18) 

The  three-term  outer  expansion  is 
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y  , 

*  *  17/  "STr  TTs'^'^TTr 
li  y  +(z-c)  /x  + 

y  ^2^^^ 

^  [x^+y^+U^^oV^ 

^71  X2<^> 

“'■  Tir  1  ,  2^  2_^,  .x2,3/2 

J-1  [x  +y  +{z-c)  1 

p/2  ^00  ^  00 

■  t7j-y/.  =“  i  L  ‘ 


X  +y  +(z-C) 


X  +y  +{z-?;) 


(4,  ^+<^-1  ^■24>i  Ht) 

'^Im  ^ly  ^In  ^Imy  1 


■  <x-m)cos  e-t-(z-n)sin  9  ^e-dm-dn 
y^+t  (x-m)cos  e+(z-n)sin  0J^ 


(7.19) 


Take  inner  variables,  the  two-term  inner  expansion  of 
the  three-term  outer  expansion  in  outer  variable  is  found  to 


<t>o+4'i+<('2'^  X 


^1^^^  -1  y  ^2^^^  -1 

-  — = - i-3n  JL  _  _4 - fan 


^  -  4— tan-"  2: 

X  2Tr  X 


^  yy2('^)  I  x‘X^(z) 

'*■  Zir  2~2  '*■  W  T~2~ 
X  +y  X  +y 

/V2  r“  r~ 

sec  6/  / 

■  TT/j  J-a>  j-00 


•  — =■  I  sec  6/ 
4ir^J-TT/2 

d6*dm*dn 

-mcos  e+(z-n)sin  6 


(7.20) 


Next,  consider  the  inner  solution.  After  a  long  and 
tedious  algebraic  manipulation,  the  mapping  inversion  formula 
gives  Z  in  terms  of  large  C  to  three-term: 
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where 


Next, 


I+Cq-  2  . 

'.U0+A1+A2-A0  j  'c,  +  ^^^Osin  a  a 

,2sin  a 


+  A  y-2— +  {A^jI-CQ-bQ-(bQCiCOt  a 

2sin  a 

+  bj^sin  Oi)sin  a  -  (cQbj^+Cj^bQ)CQl  -  Aj^Cc^+b^) 

-  i  [A^sin  a-  A^c^  cot  o+bj^sin  a) 

+  A^sin  a{cQbj^+Cj^bQ)  +  A^^sin  aJ)log  c 

+  4  AQ(Mj^4-iRj^)  +  A^CM^-iRQ)] 

+  ~  {-c^b^d+c^bQ)  +  T^sin^  a  -  tsin  al(l4-Cnb„)b 


”0  0 

♦  CoT„l)+0(^) 


0‘'0'“0 


(7.21) 


V2 


4-  b„^/-l 


- ^ —  +  On'*'n 


6.  b|4-c. 

''l  *  ^‘5;  -  E^> 

5, 

bn 

Ml  =  =1  -  d+Co^o’  ■  '=O*’l*'=l'=0”^‘’«  "’0*^' 

^nb/jC. 

Ri  =  "  b,sin  a)  log  2  4-  b,sin  a 

1  sin  a  1  1 


/-2 - 


(7.22) 


express  r,  in  terms  of  large  Z  to  three-term: 


C,  A  c.(l+c^  )  - 

!C'  %°.-i  °  <^i1) 

0  2sin  a 


*0  *0 


« 

-  { (e  -  Cq  (Cghj^+Cj^bQ)  -  (b^Cj^cot  a+bj^sin  a)  sin  a 

“  {CQ+bg)  -  i  IcQbj^+Cj^bQ)sin  a  -  (b^Cj^cot  a 

+  bjSin  a)c„  *  ^  sin  -I  '  +  ii^)  (e-"“-bo)  > 

•(log  Z  -  log  Ajj)  -  e“^®(MQ-iRQ) 

"  --5I^>  '•■^“-'>o>  -  ^*51^'  ‘V^'o’ 

A  A 

+  +  -^  Mq  +  i  (Ri'a^  ^ 

+  ^(e“^®-bQ)^*(log  Z-log  A^) 

+  ^  te“"“(Mo-iRo){e'"«-bo)  +  Cobod+c^b^) 

-  TgSin^  ot  +  iCgT^sin  a+(l+CQbQ)bQl  +  Q{^)  ^3) 

z 


If  equation  (7.23)  is  used,  the  second-order  inner 
velocity  potential  expanded  with  respect  to  large  Z  to  three- 
term  can  be  obtained.  Take  outer  variables,  one  obtains: 


Cl 

e(4>Q+4‘]^)  -v  X  +  y  gYn~a  ”  E(AQ+Aj^)sin  a*tan"-^-J 


Ci^  ycgd+c^j^) 

—T”  X  +  - . .  "■  +  eA-[(b«c,cot  a 
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c,  . 

+  b^sin  a)CQ  -  (CQb^+CjbQ)sin  a+^jjp^(Cj,+bQ)  ] 

^0 

+  [  Cq  a+bj^sin  a)  sin  a  -bj^l 


2*  2  . 

2  2.  ^  ^0  “ 


•log  (x^+y^) - 2-_ — log  (x^+y^)  +  ytan  ^  ^1 

x'^+y'"  ^  * 


•sin  a  +  (xtan”^  x  "  2  (x^+y^)l tCQ+b^) } 

^2.  2 

e  A  X 

“  — n  a-RQCQ)sin  a  -  (l+CQbQ)CQbQ 

X  +y 

2  2 
+  T^sin  a  +bQ(bQ+CQ)  -  sin  a ‘log  eAq] 

c^A  \ 

+  — (  KoCj+RjSin  o-  CjTj  -  (l+Cjbjlbj 


-  (bQ+CQ)  log  eAQjsin  a 


(7.24) 


This  is  the  three-term  outer  expansion  of  the  two-term  inner 
expansion. 

If  equations  (7.18)/  (7.20)  and  (7.24)  are  used/  the 
strength  of  singularities  ^2'  ^2  obtained. 


VII.  5.  Third-Order  Free-Surface  Expression 
The  equation  (7.3)  gives: 


^'^I'^lyy  ^2y  "  '^Ix’^lx 


(7.25) 


where  the  radiation  condition  gives  the  lower  limit  of  the 
integral  to  be  -  «. 


f 
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The  equations  (4.11)  and  (4.19)  give  (Appendix  B) 


<J)j^(x,  ±0,  z)  =  ±|-  Y]^(z) 


in  the  wake 

on  y=±0  (7.26a) 


This  function  is  well  defined  and  can  be  differentiated  with 
respect  to  x  and  z  directly  at  y=0  to  give; 


'l»lx(x»  ±0/  z)  =  0 


x%0 


(7.26b) 


♦iz'*'  <  0 


and  frcHQ  Appendix  D,  it  gives 


inside  the  wake 
outside  the  wake 


(7.26c) 


4* 


ly 


(X  +0  z)  -  ^  f  ^  r  1  . 

(X,  -0,  z)  -  4Ttf  dc  ^z-c  xiz-e) 

J  (7.2 


(7.26d) 


Substitute  equation  (7.26)  into  (7.25)  to  give; 


♦lyy  +  *2y> 


on  y=0,  x<0  (7.27a) 


Only  this  upstream  condition  is  needed.  However,  if  equations 
(7.26a)  -  (7.26c)  are  used,  one  obtains; 


‘^lyy  '^Ixx  ^Izz  ^ 


for  y=0,  x<0  (7.27b) 


Finally,  second-order  free-surface  deflection  becomes; 
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■/: 


«2y-dx 


on  y=0,  x<0 


(7.27c) 


So  the  three-tern  outer  solution  is: 


Ho+n 


1  f  1  2  2 

-|z-c|-logj 


♦2y-<Jx 


on  y*0,  x<0 


(7.28) 


Now,  consider  the  three-term  outer  expansion  of  the  two- 
term  inner  expansion  for  the  free-surface  deflection.  Set 
in  equation  (7.21a).  Then  separate  the  real  and  imaginary 
parts  and  eliminate  the  free-surface  elevation  for  large 
C  is  thus  obtained. 

In  order  to  perform  the  free-surface  matching,  equation 
(7.27c)  must  be  evaluated  near  the  line  of  singularities. 
Unfortunately,  due  to  the  mathematical  difficulty,  this  inner 
limit  has  not  yet  been  obtained.  However,  if  this  limit  is 
evaluated,  then  its  matching  with  the  three-term  outer 
expansion  of  the  two- term  inner  expansion  for  free-surface 
will  provide  an  equation  to  relate  the  depth  of  submergence 
and  second-order  parameters  a^^,  bj^,  c^^,  and  thus  give  an 
unique  second-order  solution. 


VIII.  WEDGE  SHAPE  PLANING  SURFACE 


VIII.  1.  Introducation 

In  the  previous  chapters,  it  is  assumed  that  there  is 
no  deadrise  angle  and  the  planing  surface  takes  the  simplest 
form  -  a  flat  plate. 

In  the  present  chapter,  the  theory  will  be  extended  to 
include  the  case  where  the  deadrise  angle  is  different  from 
zero.  This  will  be  the  case  for  a  V-shape  planing  surface. 
It  is  also  shown  that  the  theory  can  also  be  carried  over 
to  any  kind  of  planing  shape  if  there  is  no  camber  in  the 
longitudinal  direction, 

VIII.  2.  Formulation  of  the  Problem 

Let  6  be  the  deadrise  angle.  The  trailing  edge  is 
lying  on  a  plane  normal  to  the  incident  flow. 

For  an  ideal  irrotational  flow,  there  exists  a 
velocity  potential  (j)(x,  y,  z)  which  satisfies  the  continuity 
condition: 

(L)  V^c})  =  0  in  fluid  (8.1) 

The  Bernoulli  equation  on  the  free  surface : 

(FI)  +  4)^^  =  1  on  (8.2) 
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and  the  kinematic  free-surface  condition: 

X  X  y  z  1 


(8.3) 


Next  consider  the  body  condition.  Define  the  wedge 
surface  by  the  equation 


y  =  S(x,  z) 

or  S  =  S(x,  z)  -  y  =  xtan  o  +  ztan  6  -  y  =  0 


(S.4a) 

{8.4b) 


Then  the  body  condition  is: 


M=o  =  d!S  +  <^S  +  <j)S„ 

3n  ^x  X  y  y  22 


(8.5) 


Substitute  (8.4b)  into  (8.5)  to  give: 


(B)  <(>jjtan  a  -  (})y  +  <{i^tan  6  =  0  on  body  (8.6) 


The  radiation  condition  at  upstream  infinity  gives: 


(R) 


♦x  * 

♦y  '  »; 


as  x-^-“>,y->— ' 


(8.7) 


The  boundary-value  problem  for  the  wedge  shape  planing 
surface  is  thus  mathematically  formulated  by  (8.1),  (8.2), 
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(8.3),  (8.6),  and  (8.7) . 

Let  the  aspect-ratio  AR  be  defined  as 


AR  = 


(Span  Projected  on  tlie  Undisturbed  Free  Surface) 

Wetted  Planing  Area 

(8.8) 


VII.  3.  Par-Field  Boundary-Value  Problem  and  Solution 
Far  away  frora  the  planing  surface,  the  detail  of  the 

body  shape  is  lost  and  the  work  done  in  the  previous 
chapter  for  the  planing  plate  is  good  for  the  wedge-shape 
planing  surface  as  well  with  only  one  modification,  namely, 
h  is  not  constant  any  more  but  a  function  of  z. 

VIII.  4.  Near-Field  Boundary-Value  Problem 
The  location  of  trailing  edge  is  given  by 


y  =  h(z)  =  hg  +  ztan  9  (8.9) 

where  h^  is  the  height  of  the  trailing  edge  at  the  middle 
section  above  the  undisturbed  free  surface. 

Let  the  inner  variables  be 


X 


x4-h(z)cot  g 
e  ' 

y-h (z)  . 
e  ' 


Y 


Z  =  z 


(8.10) 
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The  inner  boundary- value  problem  is  t 


(L) 

♦xx  * 

‘J'yy 

^  =  0 

in  fluid 

(8.11) 

(FI) 

♦  ^  + 

Y 

.  e\^  =  1 

on  S^ 

(8.12) 

(F2) 

r*x“x 

- 

Y 

.  EtjNj  =  0 

on  S^ 

(8.13) 

(B) 

»xtan 

a  - 

$Y  +  6 

=  0  on  body 

(8.14) 

Assume  the  existence  of  asymptotic  series  for  and  N  as 
given  in  (2.22). 

A.  First-Order  Inner  Solution 

The  first-order  boundary  value  problem  is  exactly  the 
same  as  given  in  section  III.  1.  Therefore,  the  solution  is 
given  by  (3.10)  and  (3.19).  Now  the  quantity  h  is  a  function 
of  z. 

B.  Second-Order  Inner  Solution 
The  boundary-value  problem  is 


(L) 

ill 

^IXX 

^IYy  =  0 

in  fluid 

(8.15) 

(FI) 

d)  A 

^OX^IX 

'^oy'^’iy  ^ 

0 

on  Sj 

(8.16) 

(F2) 

^ox^ix 

“^IX^^OX  ■ 

e$lY  =  0 

on  Sj 

(8.17) 

(B) 

oc  -  $1Y  + 

=0  on  body 

(8.18) 

It  is  noted  that 
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h(z)  =  0(A{z)) 

so  h(z)  =  0(e)  as  e-*-0  (8.19) 

The  equations  (8.9)  and  (8.19)  give: 

tan  0  =  o\l)  as  e-^0  (8.20) 

Use  of  (8.20)  in  (8.14)  gives: 

^j^jjtan  a  -  =  0  on  body  (8.21) 

It  is  thus  shown  that  the  second-order  inner  boundary- 
value  problem  has  exactly  the  same  form  as  given  in  section 
IV.  3.  The  solution  is  thus  given  by  (4.31)  and  (4.37). 

VIII.  5.  General-Shape  Planing  Surface  with  Curve  in  the 
Spanwise  Direction 

Previous  sections  show  that  the  boundary-value  problem 
for  V-shape  planing  surface  has  exactly  the  same  fonn  as 
for  the  planing  plate. 

Now  consider  the  general-shape  planing  surface  given  by 
the  equation 


y  =  xtan  a  +  S(z) 


(8.22) 


106 


Equation  (8.22)  implies  that  there  is  no  camber  in  the 
longitudinal  direction. 

Next,  it  is  to  be  shown  that  the  theory  up  to  the  second- 
order  solution  presented  in  the  previous  chapter  can  be 
carried  over  directly  to  this  general  shape  planing  surface, 
with  curve  in  the  spanwise  direction.  In  order  to  show 
this,  it  is  sufficient  to  show  that  the  body  condition  up 
to  the  second-order  solution  has  the  same  form  as  given  for 
the  planing  plate. 

The  body  boundary  condition  gives: 


■^  =  0  =  4»  tan  a  -  ;{»  +  0  S  on  y=xtan  a+S(z) 

an  X  y  z  z  ,5  23, 

Let  N  denote  a  unit  vector  lying  in  a  cross-section 
plane,  z=  constant,  the  vector  being  normal  to  the  contour 
of  the  section  cut  by  the  plane.  Then  the  body  boundary 
condition  can  be  re-written 


li 

3N 


tan  a'tfi  -41^ 
/l+tan^  a 


or 


^z 

3N  A  1 

/1+tan  a 


(8.24a) 

on  y=xtan  a+S(z)  (324^) 


In  terms  of  inner  variables,  the  left-hand  side  is  0(|-), 
because  of  the  differentiation;  the  right-hand  side  is 
0(e({))  since  S  is  0(e). 
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Let  the  solution  have  the  inner  expansion 
N 

'I>(x,y,Z)  E  4>  (X,y;Z;e) 
n  n 
n=0 

where  .i=o(4>  )  as  e-^0  for  fixed  X,y,Z  (8.25) 

n+i  n 

If  (8.25)  is  substituted  into  (8.24b),  the  body  condition 
is  clearly. 


'*0 

•3^ 


0 


on  body 


(8.26) 


To  the  second-order,  it  is  noted  that  4>j^=0(l).  The  left- 
hand  side  is  0(|-)  while  right-hand  side  is  0(e).  Clearly, 


9$1 

=0  on  body 


(8.27) 


Therefore,  to  the  second-order  inner  solution,  the  body 
condition  only  depends  on  X  and  y  explicitly  at  a  constant 
z-plane. 

The  theory  developed  for  the  planing  plate  can  thus  be 
applied  directly  to  the  general  shape  high-aspect-ratio 
planing  surface  with  curve  in  the  spanwise  direction. 


IX  CONCLUSION 


A  hitjh-aspect-raiio  planing  suzface  is  investigated 
by  the  use  of  r'atcbcd  asymptotic  expansions.  A  three 
dimensional  efft'ct  is  shown  to  provide  a  height  rof-^tence 
for  the  planing  sutface  so  that  the  location  of  the  planing 
suriace  is  uniquely  defined  with  respect  to  tlie  undisturbed 
free  surface.  The  numerical  solution  is  presented  in 
Chapter  six. 

In  the  present  theory,  it  is  assiined  that  the  chord 
length  distribution  and  its  first  derivative  with  respect 
to  the  span  are  continuous  at  the  tip  as  elsewhere.  The 
cusp  shape  distribution  does  satisfy  those  requirements. 
However,  for  practical  applications  it  is  desired  to 
include  a  rectangular  plate  or  elliptic’-shape  distribution. 
The  last  two  shapes  obviously  violate  the  assumption  of 
continuous  first  derivative  at  the  tip.  Prom  analogies 
in  aerodynamics,  it  is  believed  that  the  ti>  effect  is 
only  a  local  phenomenon,  at  least  for  the  elliptic-shape 
chord  distribution.  Thus  the  theory  presented  here  for 
elliptic-shape  should  provide  a  valid  result  tlirough  the 
whole  span  except  in  a  very  small  neighborhood  of  the  tips. 
Observation  seems  to  confirm  this  statement  even  for  a 
rectangular  plate.  If  the  tip  effect  were  shown  to  be 

28 

significant,  a  tip  correction  as  suggested  by  Van  Dyke 
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night  be  applied. 

In  order  to  obtain  a  unique  second-order  sclution,  it 
is  necesrary  to  obtain  a  third-order  outer  solution;  matching 
to  the  neat-field  and  far-field  f jree-surface  elevations 
provides  an  equation  relating  the  depth  of  submergence  and 
the  unknown  second-order-problem  parameters.  Unfortunately, 
due  to  a  mathematical  difficulty,  the  inner  limit  of  the 
third-order  outer  free-surface  solution  is  not  obtained  in 
the  present  work.  However,  the  method  of  solution  is 
presented  in  Chapter  seven. 

The  bottom  of  the  planing  surface  is  assumed  to  be 
flat  through  the  first  seven  chapters.  In  Chapter  eight, 
this  restriction  is  releueed  to  include  a  V-shape  planing 
surface.  To  the  first  order,  the  present  theory  is 
applicable  for  arbitrary  angle  of  deadrise.  However,  in 
the  second-order  solution,  it  is  implicitly  assumed  that 
the  deadrise  angle  is  small. 

Throughout  the  whole  analysis,  it  is  assumed  that 
there  is  no  camber  in  the  longitudinal  direction  so  that 
the  classic  hodograph  method  is  applicable  to  obtain  the 
inner  solution.  This  restriction  can  also  be  relaxed  in 
the  manner  as  given  by  Wu.^^ 

Gravity  is  neglected  in  the  present  theory  by  assuming 
a  high  speed  planing  surface  with  no  wave  behind  it.  In 
case  the  resistance  of  the  planing  surface  is  sought. 
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energy  which  is  carried  away  by  waves  must  be  taken  into 

account.  To  do  this,  Froude  number  must  be  considered. 

There  will  be  two  small  parameters;  one  is  the  inverse  of 

large  aspect-ratio  and  the  other  is  the  inverse  of  large 

Froude  number.  To  handle  tvfo  parameters  at  the  same  time, 

33 

the  method  developed  by  Ogilvie  may  be  applicable  so 
that  the  gravity  effect  can  be  included. 
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APPENDIX  A 


RELATIONSHIP  BETWEEN  Z-And  ^-PLANES 

A.  1.  To  the  First-Order  Parameters 

Equation  (3.17)  gives  an  expression  between  the  physical 
Z-plane  and  the  ccxnplex  C“plane  for  any  value  of  To 
perform  the  matching/  only  the  outer  limit  is  required  for 
the  inner  solution.  It  is  shown  in  section  III.  3.  that  Z 
and  C  correspond  to  within  a  real  scale  factor  far  away  from 
the  planing  surface.  The  right  hand  side  of  equation  (3.17) 
can  thus  be  expanded  in  terms  of  large  C  to  give: 


Z  'v  AQe"^“{(-.CQ+  /i-Cq^)?  +  logUl+CobQ-ibp/l-CQ^) 
-i/l-CQ^(bQlog2  +  vt^‘log(bQ-^^^))  -  c^ 


-  (l+CQbQ)«log  (l+bg)}  +  0(i) 


(A-1) 


Next,  express  C  in  terras  of  large  Z. 


ta 


Ao(”Co+i»-‘  ''0 


./1-C„  ) 


Z  - 


(-Cn+t/l-Cf,  ) 


log  C  -i/l-Cp^  tbplog  2  +  v4>q^-1  log  (bQ-«4>Q^-l)  ] 


-  Cq  -  (l+Cgbp)  log  (l+bp)}  +  0(i) 


(A-2) 


where  A,.  =  _  ^  r 

0  If  (bQ+CpT 


Now,  take  logarithm  on  both  sides  of  equation  (A-2)  and 
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note  that  log  |  f:  { =o  ( 1  ^  |  > .  Equation  (A-2)  becomes : 

ia 

C  'v.  - S - ?  '  ’7  ^  *  0<log  Z) 

Ao(-Co+i/l-Co  ) 

If  (A-3)  is  applied,  first-order  velocity  potential 
matching  gives: 

Cq  =  -  cos  a  (A-4) 

Substitute  (A-4)  into  (A-1)  to  give; 

Z  a  AqC  +  AQ(e“^“  -b^)  log  i; 

+  AQe"^“(MQ  -  iRq)  +  o(i)  for  large  C  (A-5) 

where 

Rq  =  bg/l-CQ^  -log  2  +  ^Q2-l*/l-CQ^‘log  (bp-^^l) 

(A-6) 

Express  ^  in  terms  of  Z,  one  obtains  the  two-term  outer 
limit  of  one-term  inner  solution  for  C  in  terms  of  Z: 

^  ^  -bo)*iog  -  <^0  ■  ^^0^^  + 

(A-7) 

where  Aq  is  real  and  positive. 

A.  2.  To  the  Second-Order  Parameters 
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Equation  (4.37)  gives  an  expression  between  the  Z-plane 
and  the  c-plane  for  any  value  of  Following  the  same 
procedure  given  in  Section  A.l,  after  a  long  and  tedious 
algebraic  manipulation,  one  obtains  the  two-term  outer  limit 
of  the  two-term  inner  solution  for  Z  in  terms  of  c  in  equation 
(4.39a},  and  ^  in  terms  of  Z  in  equation  (4.39c);  and  the 
three-term  outer  limit  of  the  two-term  inner  solution  for  Z 
in  terms  of  ?  in  equation  (7 . 21) ,  and  ^  in  terms  of  Z  to  three- 
term  in  equation  (7.23). 


I 


APPENDIX  B 


To  show  that  in  the  wake  region,  the  density  of  doublets 
Y(x,z)  is  a  function  of  z  only. 

In  the  wake,  the  velocity  potential  <{»  is 


_ y*Y(^>C) _ _ 


(B-1) 


Let  S  represents  the  area  of  this  integration.  Next , consider 
this  surface  integral  into  two  parts.  Let 

S  =  S*  +  S'  (B-2) 


4r 

Where  S  represents  all  the  surface  area  in  (B-1)  except  at 
a  small  circle  of  radius  r  with  center  at  (x,0,z)  and  S' 
represents  the  area  in  this  small  circle. 


4) 


y*Y(^fC) 

[(x-o^+y^+(2-c) 


2j3/2 


(B-3) 


* 

Consider  the  first  integration  S  ,  since  there  is  no  singu¬ 
larities  in  the  intergrand,  this  integral  approaches  zero 
as  y  approaches  zero.  Next,  consider  the  integral  in  the 
small  circle.  When  the  radius  r  approaches  zero,  the  area 
shrinks  to  a  point.  The  density  y  inside  this  small  circle 
can  be  considered  to  be  constant. 
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lim  y(x,z) 

r+0 


JJ 


({x-C)^+y^+(z-c)^]^^^ 


The  surface  integral  is  shown  to  be  zero  if  the  area  of 
small  circle  S'  is  excluded.  Therefore,  the  surface  of 
integration  (B-4)  can  be  extended  to  infinity  to  give; 


<|»  'V  Y  (x,z) 


[(x-5) 


y  dg  dc 
^+y^+(z-C) 


2p72 


Integration  gives : 

(|>  'V.  ~  y(x,z)  sgn  y 

Equation  (4.2')  gives 


<t>^(x,y,z)  =  0 

From  equations  (B-5)  and  (B-6> 


as  y-^ 


±0 


( 


on  y=o 


one  obtains : 


( 


Y„(x,z)  =  0 


on  y=0 


APPENDIX  C 


To  (^tain  the  inner  limit  of  outer  solution 
The  two-term  outer  solution  is  given  in  equation  (4*11). 
In  order  to  obtain  the  inner  limit  for  this  outer  solution, 
this  solution  must  be  evaluated  near  the  line  of  singularities 
and  the  equation  becooies  divergent  integral.  To  avoid  this 
difficulty,  the  method  of  integration  by  parts  is  used.  Note 
that  the  strength  of  singularities  is  zero  outside  the  planing 
body. 


Yj^(z)  =  Mj^(z)  =  Xj^(2)  =0  at  z*+l,-l  (C-1) 

Then  equations  (4.11)  beccmes: 


Take  the  inner  variables  in  {C-2)  to  give: 


i}io+4ii  '''  (eX-eHcot 


a)  - 

1  eYteH 

(eX-eHcot  a)  W  d?  °  I 


4tt  (eX-eHcot  a)  d?  ~J  ^  ^ 

1  eY+eH 

+  ^  ’  2  2  2y.  (Z) 

(eX-eHcot  a)  +{eY+eH) 
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+ 


L-. 

4w 


cX-eHcot  a 

T 


(eX-eHcot  a)‘+(£Y+eH) 


2Xj^(Z)  +  o(€Yj) 


If  equation  (4.10)  is  used,  one  obtains: 


*0**1  <eX-EHcot  o)  -  J^(eY+eH) 


Yi(Z) 

2v 


tan  ^( 


eY+eH 


'i> 


.  11- 
z-c 


eX-eHcot  o 


) 


(eY+£H)Mj(Z)  +  (eX-eHcot  a)Xj^(Z) 
(eX-eHcot  a) (eY+eH) ^ 


+  o(cX^) 


(C-3) 


'•s 


APPEie)IX  D 


To  show  that  equation  (5.16)  gives: 


1  1 


Let  S  =S-fW  represents  the  area  of  the  planing  surface 
and  the  wake.  The  second-order  velocity  potential  is 


1  rr  3nri<5»5) 


dC-dc 


((x-C)^  + 


(D-1) 


(D-2) 


Then  equation  (D-1)  becomes: 


il(x,y,2)  ■=  -  ^  //  Xi(C,c)5|-(5)dC' 

JJs+w 


Next,  differentiate  (|>j^  with  respect  to  y. 


♦ly<>= 


(D-3) 


Next,  apply  the  method  of  integration  by  parts  in  equation 
(D-3)  to  give; 
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Therefore, 


(D-5) 


123 


Substitute  equation  (0^5)  inlx>  equation  {D-4)  to  give: 


v*,„(x,y,z) 


6(0) 

J-1  '  Jo  ^  l(x-e)^+y^lR 

+  (5zi>  +r 

y‘^+{z-C)^  Jo  [y^+(z-C)^lR 

1  r  ^**^1  r  X  (z-C)  tx^+2v^+  (z-C) 


+  - jldc 

y%(z-j;) 


(D~6) 


Now,  the  integrand  in  equation  (D-6)  is  well  defined,  so 


(X  0  z)  -  -  4-  r  +  J:_] 

(x,0,z)  -  fjr/^  x(z-?)  *  z-5^ 


(D-7) 


where  Cauchy  Principal  Value  is  taken. 
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i»  aaitnacT  ^  high-aspect-ratl o  planing  surface  gliding  on  a  stream  of  an  infinitely  deep, 
incompressible,  inviscid  and  gravity-free  fluid  is  treated.  This  complicated  problem  is 
decomposed  into  two  relatively  simpler  boundary-value  problems:  1)  The  near- field 
boundary-value  problem  is  valid  only  in  the  neighborhood  of  the  planing  surface.  The 
problem  is  solved  by  the  classical  hodograph  method.  The  second-order  inner  problem  is 
also  shown  to  be  a  plane,  irrotational  flow  and  the  solution  is  obtained  by  following  | 
the  same  procedure  as  in  the  first-order  inner  solution.  2)  The  far-field  boundary- 
value  problem  is  valid  only  far  away  from  the  planing  surface.  The  first-order  outer 
solution  is  shown  to  be  a  trivial  uniform  flow.  The  outer  velocity  potential  is  defined 
in  the  whole  space  by  harmonic  continuation.  The  second-order  solution  is  then  shown  to 
be  similar  to  a  lifting-line  solution. 

The  unknown  strength  of  singularities  is  obtained  by  matching  of  the  velocity 
potential.  Then  a  matching  of  the  free-surface  deflection  provides  a  height  reference  j 
for  the  planing  surface.  The  location  of  the  planing  surface  with  respect  to  the  undis¬ 
turbed  free  surface  is  uniquely  defined.  In  order  to  obtain  a  unique  second-order  solu¬ 
tion,  it  is  necessary  to  solve  the  third-order  outer  solution.  The  detail  of  this  solu¬ 
tion  is  presented. 

A  numerical  solution  for  a  planing  plate  of  arbitrary  angle  of  attack  is  pre¬ 
sented.  A  downwash  correction  is  also  included.  It  is  shown  mathematically  that  the 
present  theory  can  be  applied  to  V-shape  or  general-shape  planing  surfaces  with  curva¬ 
ture  in  the  spanwise  direction. 
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